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ABSTRACT 

A number of radio interferometers are currently being planned or constructed to observe 21 cm 
emission from reionization. Not only will such measurements provide a detailed view of that epoch, 
but, since the 21 cm emission also traces the distribution of matter in the Universe, this signal can be 
used to constrain cosmological parameters at 6 < z < 20. The sensitivity of an interferometer to the 
cosmological information in the signal may depend on how precisely the angular dependence of the 21 
cm 3-D power spectrum can be measured. Utilizing an analytic model for reionization, we quantify all 
the effects that break the spherical symmetry of the 3-D 21 cm power spectrum and produce physically 
motivated predictions for this power spectrum. We find that upcoming observatories will be sensitive 
to the 21 cm signal over a wide range of scales, from larger than 100 to as small as 1 comoving Mpc. 
Next, we consider three methods to measure cosmological parameters from the signal: (1) direct fitting 
of the density power spectrum to the signal (this method can only be applied when density fluctuations 
dominate the 21 cm fluctuations), (2) using only the velocity field fluctuations in the signal, (3) looking 
at the signal at large enough scales such that all fluctuations trace the density field. With the foremost 
method, the first generation of 21 cm observations should moderately improve existing constraints on 
cosmological parameters for certain low-redshift reionization scenarios, and a two year observation 
with the second generation interferometer MWA5000 in combination with the CMB telescope Planck 
can improve constraints on fl w (to ±0.017, a 1.7 times smaller uncertainty than from Planck alone), 
fl m h 2 (±0.0009, 2.5 times), h 2 (±0.00012, 1.5 times), (±0.003, 3 times), n s (±0.0033, 1.4 
times), and a s (±0.003, 2.7 times). Larger interferometers, such as SKA, have the potential to do 
even better. If the Universe is substantially ionized by z ~ 12 or if spin temperature fluctuations are 
important, we show that it will be difficult to place competitive constraints on cosmological parameters 
from the 21 cm signal with any of the considered methods. 

Subject headings: cosmology: theory - intergalactic medium - reionization 


1. INTRODUCTION 

The reionization of the Universe involves many poorly 
understood astrophysical phenomena such as the forma¬ 
tion of stars, the escape of ionizing photons from star¬ 
forming regions, and the evolving dumpiness of the gas 
in the intergalactic medium (IGM). However, reioniza¬ 
tion imprints signatures onto 21 cm emission from high- 
redsliift neutral hydrogen, as will be studied with the in¬ 
struments PAST, LOFAR, and MWA, in a manner that 
is sensitive to these processes. 4 Moreover, the 21 cm 
emission encodes information pertaining to fundamental 
cosmological parameters. Due to all the overlying astro¬ 
physics, it is uncertain whether or not 21 cm observations 
can be competitive with other cosmological probes. 

Several authors have discussed using the 21 cm signal 
from reionization to study cosmology, in add ition to map¬ 
ping o ut the epoch of rei o nization (EOR1: llScott fc Reesl 
1990j_|Tozz^^nj2000l llliev et all 120021) . Recently, 
RarkanT'fcT!oeb|l|2005ari show that redshift-space dis¬ 


tortions from peculiar velocities allow for the decompo- 
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sition of the observed 21 cm 3-D power spectrum into 
terms that are proportional to /i°, /i 2 and /i 4 , where ^ is 
the cosine of the angle between a mode k and the line- 
of-sight (LOS). In principle, this decomposition makes 
it possible to separate the contribution from reionized 
bubbles from that owing to a fundamental cosmological 
quantity, the linear-theory density power spectrum. 

Even if the signal from t he ioniz ed bubbles dominates 
over the cosmological one, iNusse 3 120051) shows that one 
can look for an asymmetry between the 21 cm signal 
measured in depth and that measured in angle. The pres¬ 
ence of this asymmetry may imply that the cosmology as¬ 
sumed in the analysis is incorrect [the Alcock-Paczynski 
(AP) effect]. This effect could help furthe r constr ain Q.„, 
and h, as well as dark energy models (N usseri [2005 1. 
It might be possible to distinguish the AP effect be¬ 
cause it creates a /i 6 dependence in the 3-D power spec¬ 
trum, which is distinct fro m the behavior that arises from 
velocity-field effects alone Jn usse iHM! lBaxka/nalll2005l) . 

For both the /it decomposition of the 21 cm power spec¬ 
trum and the AP effect, the feasibility of inferring cosmo¬ 
logical parameters using future surveys depends on how 
sensitive these surveys are to deviations from sp heric al 
symmetry in the 3-D power spectrum. iMoraled ( 2005 1 
suggests that 21 cm observations should spherically av¬ 
erage fc-modes over a shell in Fourier space to increase the 
signal-to-noise ratio. In addition to losing the angular in¬ 
formation contained in the signal, such averaging would 
significantly bias upcoming measurements: The power 
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spectrum is not close to spherically symmetric and 21 
cm interferometers will be most sensitive to modes ori¬ 
ented along the LOS. Array design also factors into the 
sensitivity to certain features in the signal. The first gen¬ 
eration of EOR arrays are still being planned, and so it 
is important to understand different design trade-offs. 

Authors have considered the 21 cm emission sig¬ 
nal in several limits, such as when the spin temper¬ 
ature fluctuatio ns are still important and the HII re¬ 
gions are small (iBarkana fc Loebl l2005bl iShaoiro et alJ 
B or when the spin temperature is much larger 
than the cosmic microwave background (CMB) tem¬ 
perature (iZaldar riaga et aHl2004 iFurlanetto et al'KQoi 
iSantos et alJ 1200511 . In reality, we do not know how 
quickly X-rays from the first stars, black holes and shocks 
will heat the gas and how quickly the spin temperature 
can couple t o the gas temperature via the W outhuysen- 
Field effect flWouthuvsenl 11952: iFieldl T958 L Previous 
work suggests that t hese proce s ses act quickly after the 
fir st stars turn on (lohlEooil j Venkatesaimff^t 1 ] lf200l 
lOhen fc Miralda-Escudel 12004 lOiardi fc MadaiJ 200311 
Moreover, it is expected that quasars or stellar sources 
contribute the bulk of the ionizing radiation. If this is the 
case, reionization will be a patchy process, in which the 
HII regions grow around these sources as the ionization 
fraction increases. 

Upcoming observations will be most sensit ive to lo wer 
redshi fts (z ~ 6—12) during reionization (Bo wman et al.l 
UMa). At these low redshifts, it is likely that the spin 
temperature is greater than the CMB temperature and 
that the ionized fraction is of order unity and very patchy. 
This is the regime that we consider for much of this pa¬ 
per. It is also possible that the ionization fraction is near 
zero for a period at these low redshifts, which will facil¬ 
itate cosmological parameter estimation. We consider 
this case as well. 

The organization of this paper is as follows. In we 
make physically motivated predictions for the form of 
the 3-D 21 cm power spectrum. We th en ge neralize the 
detector noise calculation of iMoralel 1 2005 1 to a power 
spectrum that is not spherically symmetric m and in¬ 
corporate realistic foregrounds into this sensitivity cal¬ 
culation (Q. These calculations allow us to estimate 
the sensitivity of upcoming interferometers to the 21 cm 
power spectrum (0. We conclude with a discussion of 
how the 21 cm signal can be used to measure fundamental 
cosmological parameters as well as a Fisher matrix anal¬ 
ysis to estimate how precisely future observations can 
constrain these parameters ( 99 - 

In our calculations, we assume a cosmology with Q m = 
0.3, U A = 0.7, S2 b = 0.046, H = 100/ikms -1 Mpc -1 
(with h = 0.7), n = 1, and cjr = 0 .9. consistent with 
the most recent determinations dSoergel et al2003) . All 
distances are measured in comoving coordinates. 


2. VELOCITY FIELD EFFECTS 

The difference between the observed 21 cm brightness 
temperature at the o bserved freq uency v and the CMB 
temperature today is (IFieldl1195H ) 


T b (x) 


3c 2 h P A 10 n H {x) a 3 [ T s (x ) - T CM b(-z)] 
32nk B T s (x)is 0 


dr 


dv ’ 
(1) 


where c is the speed of light, kp is the Boltzmann con¬ 
stant, hp is Planck’s constant, a = 1/(1 + z), x is the 
spatial location, Tomb is the the CMB temperature, 
Aio = 2.85 x 10 -15 s -1 is the spontaneous 21 cm tran¬ 
sition rate, T s is the spin temperature, uq = 1420 MHz, 
and riH is the number density of neutral hydrogen. The 
factor \dr/dv\ accounts for the Hubble flow as well as 
peculiar velocities. If we take the average of equation m 
we find 


T b ss 26 x H 


T s — Tcmb \ ( H b /i 2 


T s J V0.022 J \V, m h 2 10 


0.15 l + z\ 1/2 

■ I mK, 


( 2 ) 

where x B is the global neutral hydrogen fraction. Fluc¬ 
tuations around T b are at the tens of mK level on mega¬ 
parsec scales. 

To calcu late dr/dv , we relate comoving distance to 
frequency llBharadwai fc AlfeOOd 8 ) 


r = 


"0( l-v r /c) 


eda 
2 H(a) 


( 3 ) 


where v r is the LOS peculiar velocity and H(a) is the 
Hubble constant. Differentiating this expression, we find 


dr c i 1 dv r 

dv o?vqH Ha dr ’ 


( 4 ) 


where we have dropped terms of order [(77a) -1 dv r /dr ] 2 
and v r /c. In the limit T s Tcmb, fluctuations in the 21 
cm brightness temperature at x can be expressed as 


^§^ = ( 1-ii [! + $*(*)]) (! + *(*)) 
J 6 


X 


L _1 dv r (x) \ 

V Ha dr ) 


- Xh, 


( 5 ) 


where Xi = 1 — x B is the global ionized fraction, 6 X is the 
overdensity in the ionized fraction and 5 is the dark mat¬ 
ter overdensity (at the scales and redshifts of interest, the 
baryons trace the dark matter). We define the normal¬ 
ized temperature as T b = T b /x B ■ In Fourier space, since 
the linear theory velocity at redshifts where dark energy 
is unimportant is v(k,z) = —iHak5i,/k 2 , the peculiar 
velocity term is 5 V = ( Ha ) -1 dv r /dr = —h 2 5l where 
fi = k ■ n, the cosine of the angle between the wavevector 
and the LOS, and L denotes the linear theory value. 5 If 
we keep terms to second order in {5, S L } } the brightness 
temperature power spectrum is 


T b 2 P/\r{k) = \x 2 h Pgs + Pxx - 2xh P x s + PxSxs] 

+ 2/ii 2 [x 2 h P Sl s - x H Px 6 l \ + ix 4 [x 2 h Ps l s l \ 

+ [2P X 6 6vx + P x <5„ 5„ x\ j (6) 

and we note that P xx = x 2 Ps n .s x and P x s = XiP bx s- 
In our calculations, we drop the connected part and set 
P x s x s = P x $ + P xx Pss ■ In equation El we have decom¬ 
posed the power spectrum into powers of p,; the last 


5 The velocity field at z ~ 10 is in the linear regime for k < 
5 Mpc -1 . See IWang fc Hi J <200511 for a discussion of the effect 
of the non-linear velocity field on the 21 cm signal. Upcoming 
interferometers are most sensitive to scales at which the velocity 
field is linear. 
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bracket in this decomposition has a non-trivial depen¬ 
dence on p. For notational convenience, we refer to 
the k-dependent coefficients in equation © as o, 2 
and P u 4 and the terms in the last bracket as Pf^ k )- 
IBarkana fc Loebl ll2005aft argue that the above decom¬ 
position should allow one to extract the “physics” - 
{PSlSl )- from the “astrophysics”- ( P xx and P x $). The 
terms in the last bracket in equation © were omitted 
in their analysis, but must be included if reionization is 
patchy because S x ~ 1 on scales at or below the bub¬ 
ble size. If we drop the connected fourth moments, the 
terms are given by 6 


Px 5 S v x{k) 


d 3 k' 

(2tt) 3 


(n-fc ') 2 [PxS L (k 1 ) P xS (\k-k'\) 


+ Ps L s{k') P xx (\k-k'\)], 
r AS I.' 

Px SvSv x{k) = J ^[{n-k'fPs LSL {k')P xx {\k-k'\) 


+( ‘ h k ' f { h '[w=h 

x P x s L (k') Px5 L (\k-k'\l 


( 7 ) 


where k! = (k! ■ k , ) 1 / 2 . These terms can contaminate the 
p decomposition. On large scales, Pf(^,k) does not de¬ 
pend on p and therefore will contaminate only measure¬ 
ments of P^a. As we go to progressively smaller scales, 
P/(/j,fc) contributes power to higher order terms in p. 

The evolution of the ionized fraction over a mode can 
also affect the spherical symmetry of Patj since time 
is changing in the LOS direction but not in the angular 
directions. The magnitude of this effect depends strongly 
on the morphology of reionization and is discussed in 
Appendix IaI 


2 . 1 . Model 

To model reionizaton by stellar sources, simulations 
must resolve halos at least as small as the HI cooling 
mass (~ 10 s M 0 ) as well as have boxes large enough 
to sample the size distribution of HII r egions, which can 
each reach larger scales than 10 Mpc llFnr1a,netto et a, 1.1 
120041 IBarkana fc Loebll2005bD . Recent simulations have 
made significant strides towards a ttaining these goals 
lllliev et al.l2005HKohler et al . 2005lh For the time being, 
semi-analytic models of this epoch provide a convenient 
approach for modeling reionization on the large scales 
that are relevant for upcoming observations. In this pa¬ 
per, we employ the physically motivated semi-analytic 
model described in Furlanetto et al. (2004a; hereafter 
FZH04) to calculate Pat- 

Re c ent numerical simulatio ns (e.g., fSokasia n et alJ 
12003 . 200 1 : ICiardi et a,1.1 (200l show that reionization 
proceeds “inside-out” from high density clusters of 
sources to voids, at least when the sources resemble star¬ 
forming galaxies (e.g., Springel & Hernquist 2003; Hern- 
quist & Springel 2003). We therefore associate HII re¬ 
gions with large-scale overdensities. We assume that a 
galaxy of mass m ga \ can ionize a mass ( m ga i, where ( is a 
constant that depends on: the efficiency of ionizing pho¬ 
ton production, the escape fraction, the star formation 
efficiency, and the number of recombinations. Values of 

6 These equations are exact if fluctuations are Gaussian. 


C < 10 - 40 are reasonable for normal star formation, 
but very massive st ars can i ncreas e the efficiency by an 
order of magnitude l |Bromm et al . 2002 ). 

The criterion for a region to be ionized by the galax¬ 
ies contained inside it is then / co ii > C -1 , where / co ii is 
the fraction of mass bound to halos above some mini¬ 
mum mass m m i n . We assume that this minimum mass 
corresponds to a virial temperature of 10 4 K, at which 
point hydrogen line cooling becomes efficient. The func- 
tion fi- n n depends on t he assumed halo mass function. 
iFurlanetto et alJ (2005 ) find that the choice of the mass 
function has an insignificant effect on the FZH04 model. 
Here we use the Press-Schechter m ass fu nctio n. In 
the extended Press -Schechter model llBond et alJll99lt 
lLacev fc ColeifTool the collapse fraction of halos above 
the critical mass m m j n in a region of mean overdensity 
S m is 


/coll = erfc 


Sr. - S„ 


V 2 [ CT mi„ -T 2 (m,z)] 


( 8 ) 


where cr 2 (m,z) is the variance of density fluctuations on 
the scale m, a 2 nin = cr 2 (m m i n , z) and S c ~ 1 . 686 , the 
critical density for collapse. With this equation for the 
collapse fraction, we can write a condition on the mean 
overdensity within an ionized region of mass to, 

S m > S B (m, z) =5 C - V2K(C)[aP n - a 2 (m.,z)] 1/2 , (9) 

where A'(£) = erf -1 (l — C -1 )- 
FZH04 showed how to construct the mass function of 
HII regions from S B in a manner analogous to that of the 
halo mass function llPress fc Schechteifl974t iBond et all 
[m2)- The barrier in equation © is well approximated 
by a linear function in a 2 , S B ss B(m) = Bo + B\a 2 {m), 
where the redshift dependence is implicit. In tha t case. 
the mass function has an analytic expression ( Sheth 

ITflSai : 


n(m) 


[2 p 

diner 

B 0 

exp 

B 2 (m) 

V 7T TO 2 

d In to 

a(m] 

2 ct 2 (to) 


, (10) 


where p is the mean density of the Universe. Equation 
m gives the comoving number density of HII regions 
with masses in the range ( m,m + dm). The crucial 
difference between this formula and the standard Press- 
Schechter mass function occurs because S B is a decreas¬ 
ing function of m. The barrier becomes more difficult to 
cross on smaller scales, which gives the bubbles a char¬ 
acteristic size. 

We can calculate P xx and P x s for the FZH04 
model with the sem i-analytic method described in 
iMcQuinn et al.l (<2005l). It is more difficult to calculate 
P„x z . IMcQuinn et a,1 ( 2005 ) showed that it is not nec¬ 
essary to consider bubble substructure in this analytic 
model when constructing P x s', ~ only the size of a bubble 
and the average density within a bubble are important. 
Since we can ignore bubble substructure, this implies 
that P x s — > P X $ L when the effective HII bubble radius 
reaches linear scales. This happens in the FZH04 model 
when Xi > 0.25. In the opposite limit, when the bubbles 
are very small, this term is subdominant to density fluc¬ 
tuations and so has a small effect on the power spectrum. 
We set P x s L = P X .S for all times. This will result in an 
overestimate of this term when Xi is small and therefore 
an underestimate of P^. 
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For our calculations in this section, our objective is not 
to model the 21 cm power spectrum for different reioniza¬ 
tion scenarios and discuss morphological differences. In 
the cont ext of the FZH04 model , this has been done in 
FZH04, iFurlanetto et alJ l|2004blh and iFurlanetto et alJ 
n®. Instead, we restrict ourselves to one parameteri¬ 
zation of this model, setting £(z) = 12, in order to illus¬ 
trate the effect of redshift-space distortions on Pat- For 
this parameterization, the EOR spans roughly the red- 
shifts 8—15. In reality, £ will have some time dependence, 
and it may even have a very complicated evolution. For¬ 
tunately, we find that the parameterization £ = 12 is 
representative of the FZH04 model: while varying the 
function f(z) will change Xi{z), if we identify the same 
ionization fraction for different values of f, the values of 
P xx are quite similar (FZH04). 

Of course, reionization might proceed differently than 
in this analytic model. The parameter ( may depend 
on the mass of the dark matter halo. This can have 
a modest effect on the characteristic size o f the bub¬ 
bles and the large-scale bubble bias iFurlanetto et alJ 
120051) . Also, recombinations might play a larger role 
in shaping the morphology of reionization. Employing 
reasonable assumptions for the dumpiness of the IGM, 
IFurlanetto fc Ohi ( 2005 1 show that the effect of recombi¬ 
nations in the FZH04 model is only important at Xi > 
0.7, increasing in importance as Xi —► 1. The presence 
of mi ni-halos m ay make recombinations more important 
(IShaniro et al.ll2fl03l) . but recent work su ggests that this 
may not be the case (ICiardi et al .11200 511 , So long as re- 
combinations are not important . IFurlanetto et alJ (I2005D 
showed that any model where the ionizing sources trace 
the collapsed regions must have a qualitatively similar 
bubble size distribution to the FZH04 distribution. 

Figure^plots for f = 12 the components of the dimen¬ 
sionless power spectrum, k 3 Pat/^it 2 , that have different 
p dependences. The thick solid, thick dashed and thick 
dot-dashed curves indicate the P M o , P M 2 , and P M 4 terms, 
respectively. The p dependence of P/( Mj fc) is nontrivial. 
The three thin solid curves indicate the total contribu¬ 
tion from P/( Mj fc) f° r P 2 = 0.0, 0.5 and 1.0 (in order of 
increasing amplitude). 

At Xi = 0.1, P^4 is the largest of the p- terms (see 
Fig. ^ right). This is because the neutral regions are 
underdense on average, and, at small ah, this results in a 
suppression of P^o and P M 2 . The p decomposition of the 
signal can be much different at larger ionization fractions. 
For Xi = 0.7, Pp 0 is much larger than the other terms at 
most scales (see FigEl left), and P M 2 becomes negative at 
small k because the neutral gas is underdense on average. 
Finally, at smaller scales than the bubbles size, P/( M ,fc) is 
larger than P^2 and P^ 4 , and is even larger than P^o at 
k k, 3 Mpc -1 . 

The evolution of the spherical symmetry as a func¬ 
tion of Xi is non-trivial. When the ionization fraction is 
small, the redshift space distortions are important on all 
scales. In the opposite limit, when the ionization frac¬ 
tion is large, the redshift-space effects are less important 
since the bubble-bubble term P xx , which enters through 
P M o, dominates the signal (Fig. Q left). The evolution 
of the angular symmetry of the signal is illustrated in 
Figured which plots the contours of constant k 3 Pat(&) 
for Xi = 0.1 and Xi = 0.7 ( right and left panels, respec- 



Fig. 1.— The //. decomposition of the signal (see equation [7J for 
Xi = 0.1 and 0.7, corresponding to z = 13.5 and 9 in the £ = 12 
model. The thick solid, thick dashed and thick dot-dashed curves 
show P M o, P„ 2 , and 4 , respectively. The three thin solid curves 
show calculated using equation Q with p, 2 = 0.0, 0.5 and 

1.0 (in order of increasing amplitude). 


tively). When Xi = 0.7, the signal is fairly symmetric at 
smaller fc-modes than the bubble scale. At larger values 
of k or at small ionization fractions, density fluctuations 
dominate the signal, and the power spectrum can be very 
asymmetric. Because of this, it may be possible to de¬ 
termine the characteristic bubble size by observing the 
angular dependence of Pat- 

In Figure El we plot the signal [k 3 PAr(k, p)/2 tt 2 ] for 
modes with p 2 = 0.0,0.5, and 1.0 (the solid, dot.-dashed, 
and dashed curves, respectively) at four times during 
reionization. Between Xi ~ 0.0 —0.1, density fluctuations 
dominate the signal and the signal is very asymmetric. 
By Xi = 0.5, the neutral fraction fluctuations contribute 
most of the power on large scales. When these fluctu¬ 
ation dominate, the 21 cm power spectrum develops a 
“shoulder” on scales near the characteristic bubble size. 
This feature moves to larger scales as the bubbles grow. 

In $3 we show that upcoming interferometers will be 
more sensitive to modes with certain orientations rela¬ 
tive to the LOS. It turns out that arrays that are very 
compact (i.e. have most of their antennae within a 1 km 
region), such as MWA, are most sensitive to fc-modes 
that are oriented along the LOS. The fact that these 
modes have more power enhances MWA’s sensitivity. 
Conversely, it is difficult to separate the terms P^o, P M 2 
and P„4 with observations that are most sensitive to the 
modes along the LOS. 

3. SENSITIVITY TO THE 21 CM POWER SPECTRUM 

In this section, we summarize how to calculate the vari¬ 
ance on interferometric observations of the 3- D 21 cm 
power spectrum. Our calculation follows that of 'Moralesl 
( 2005 ). and extends their calculatio n to capture the an- 
gular dependence of the 3-D signal. I White et all (|1999D 
and lZa,lda,rria,ga et a,U (l2004f> do a similar interferomet¬ 
ric detector noise calculation, but for the angular power 
spectrum. For 21 cm observations, the 3-D power spec¬ 
trum is more interesting than the angular power spec¬ 
trum owing in part to the the p dependence of the sig¬ 
nal. In addition, the 3-D power spectrum will allow us 
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log[k/(.01 Mpc ^ cos(6) 


Fig. 2.— Contours of constant & 3 Pat W for the same signal as in 
Figure^ The horizontal axis is the LOS direction, and the vertical 
axis is the transverse direction. The coordinate transformation 
(/c_|_, /c| |) —* log(/c/(0.01 Mpc _1 ))(sin(0), cos(0)) preserves circles of 
constant power. 



k (Mpc ^ 

Fig. 3.— Dimensionless power spectrum A: 3 Pat/27t 2 for p? 
equal to 0.0, 0.5, and 1.0 (solid, dash-dot, and dashed lines, re¬ 
spectively) for four times during reionization in the £ = 12 model, 
corresponding to z = 9, 9.8, 13.5 and 20 (in order of decreasing Xi). 
The signal is most asymmetric at scales where density fluctuations 
dominate. 


to measure many more independent modes. 

The 21 cm signal will be observed with radio interfer¬ 
ometers, which measure visibilities. The visibility for a 
pair of antennae, quantified as a temperature, is given 

by 

V(u t v, v) = j dft Ari(n,!/)4(n)e 2 "(*) 4 , (11) 

where (it, v) is a vector that gives the number of wave¬ 
lengths at frequency v between the antennae and A v (n) 
is the contribution to the primary beam in the direction 
h. Here, we are working in the flat sky approximation; 
this is adequate since upcoming experiments are most 
sensitive to angular modes with wavelength 9 < 0.1 ra¬ 
dians. 

We assume that the visibilities are complex Gaussian 
random variables, such that the likelihood function of the 
covariance matrix Cij = {V* Vj) for n visibilities, where 


the asterisk indicates a complex conjugate, is 


C(C) 


1 

7r n detC 


exp 


E^T 

id 


( 12 ) 


Because the visibilities are complex and V (it, v, v) = 
V(—u,—v,v)*, when counting the number of indepen¬ 
dent u-v pixels we restrict ourselves to the half space. In 
this section, it will be more convenient to work with the 
Fourier transform of V in the frequency direction. This 
is just a change of basis of V and C in equation o. 

For upcoming arrays, C will be dominated by the 
detector noise on most scales. The RMS detector 
noise fluctuation per visibility of an antennae pair af- 
ter observing for a tim e to in one frequency channel is 
iiRobifs fc wiisoniliinnl 


av n 


\ 2 t 

/v sys 

A e \J Av t 0 ’ 


(13) 


where T sys (v) is the total system temperature, A e is the 
effective area of an antenna, and Av is the width of the 
frequency channel. 

For an observation with bandwidth B , where B Av , 
if we Fourier transform the observed visibilities in the 
frequency direction, we then have a 3-D map of /(it) = 
f dvV (it, v, v) exp[ 27111177 ], in which it = ui+v j+rj k and 
r/ has dimensions of time. If we perform this transform 
on just the detector noise component V N of the visibility 
map V, we have 


B/Av 

I N (u)= ^ V N {u, v, v^ exp[2nivir]] Av 
i =1 


B/Al> 

= ^2 V ,N (u,v,Vi) Av, 

2=1 


(14) 

(15) 


where we have absorbed exp[27rw?7] into a new variable 
V' N that has the same RMS as V N and the frequency 
channels tq,..., vb/Av ar e spaced Av apart. It follows 
that the d etector n oise c ovariance matrix for a single 
baseline is l Moralesil2005f l 


/ B / Av 


= ( E V ' N ( U ii V m) Av 


\ ra=1 
^ B/Au 


X 


E V ,N (uj,vi)Av 


1= 1 


= B Av (AV n ) 2 Si 


A 2 BT, 


sys 


K. 


B to 


(16) 

(17) 


To reach equation d, we note that different 
V' N (ui,v m ) are uncorrelated, and equation CE3 follows 
from equations and Note that equation CO 

only depends on B and not on Av: finer frequency reso¬ 
lution comes at no added cost. 

We now estimate the average observing time tk for an 
array of antennae to observe a mode k as a function of 
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the total observing time to (note that there is an isomor¬ 
phism u <-> k and that 2nu±/x = fcj_, where x is the 
conformal distance to the emission). At any time, the 
number density of baselines that can observe the mode 
u is We assume this density is rotationally in¬ 

variant and define 9 to be the angle between k and the 
LOS. Integrating n(u±)dudv over the half plane yields 
A/base ~ IVf nt /2, where N an t is the number of antennae. 
Since the telescope observes a region in the u-v plane 
equal to SuSv ~ A e / A 2 , each visibility is observed for a 
time 

tk ~ n i x \ k \ sin(0)/27r), (18) 


where to is the total observing time for the interferome¬ 
ter. It follows that the detector noise covariance matrix 
for an interferometer is 


C N {ki,kj) 


A 2 B T sys \ 2 Sjj 
A e ) B t ki 


(19) 


(From now on we will use k rather than u to index ele¬ 
ments in C). 

We also want an expression for the contribution to C 
owing to sample variance. For a 3-D window function 
W(ra, v) = A v {h) F^iu). if we assume that different pix¬ 
els indexed by u are uncorrelated, the covariance matrix 
of the 21 cm signal I 21 is 

C sv (k i ,k j ) = (P\u i )*I 21 (u j )) 

« % jd^u'\W{ui-u')\ 2 P^ T {u')m 


wavelength corresponding to k fits within the survey vol¬ 
ume (i.e. when 27r/fc cos((9) < y) and otherwise we set 
N c = 0. 7 

When we sum equations m and El to get C, the 
error in P]^ T (k) from a measurement in an annulus with 
N c {k : 9) pixels is 8 


5P^ T (k,e) « [C sv (k,9) + C N (k,9)] , 

(24) 

where we have defined C (k) = C(k , k). One can trivially 
derive equation (EH by calculating the Fisher matrix F 
of the P(fcj) from C (e.g. eqn. ED with C —► C), 
noting that there are N c measurements of P(k,9) and 
that 6P$r(k) = VF~ 1 k k- 

Given a model for t he data , the 1-u errors in the model 
parameters A * are (v / P ,_ 1 ii), where 


. _ y- 1 OPat 9Pat 
lj ~ (5P AT ) 2 ~d\~ ~d\~ 


(25) 


[See Appendix iBl for a useful formula for SPAr{k), the 
error in the angular-averaged power spectrum.] 

In the next section, we extend the above analysis to 
include foregrounds. The calculation in this section as¬ 
sumes Gaussian statistics, but the ionization fraction 
fluctuations on the scale of the HII bubbles will not be 
Gaussian. Numerical simulations are necessary to quan¬ 
tify the degree of non-Gaussianity introduced by patchy 
reionization. 


where we have used that (AT 21 (W)AT 21 (u)) = 

Pat{u) S 3 (u' — u) and the definition of visibility (eqn. 


Hill. We can simplify C sl further: 



P 21 (u ■) ^ B 5 

r A T\ u i) ^ °ij 

(21) 


\2 td2 

PAT(ki) Sy, 

A e x 1 y 

(22) 

where to get to equation 

ED) we pull P{ Ui ) 

out of the 


integral and use the fact that W (it) is different from zero 
in an area SuSvSr) fts A e /(X 2 B) and must integrate to 
unity within the beam, such that f d 3 u' \ W(u — u/)| 2 ~ 
(SuSvSrj)- 1 . Equation (l2H is accurate for values of |tq| 
much greater than the FWHM of W. The additional 
factor of x 2 y/B that arises in equation ED is because 
with our Fourier conventions Pat(^) = x 2 y/B Pat(u), 
where y is the conformal width of the observation. 

Over the course of an observation, a large number of 
independent Fourier cells will be observed in a region of 
real space volume Vol = x 2 y A 2 /A e . We have seen that 
the 21 cm power spectrum is not spherically symmetric, 
but is symmetric around the polar angle </>. Because of 
this symmetry, we want to sum all the Fourier cells in 
an annulus of constant (k, 9) with radial width A k and 
angular width A 9 for a statistical detection. The number 
of independent cells in such an annulus is 

N c (k, 9) = 2nk 2 sin(0) A k A 9 x (23) 

Here, (2-7r) 3 /Vol is the resolution in Fourier space. For 
our calculations, we use equation El for N c when the 


4. FOREGROUNDS 

The foregrounds at A = 21 (1 + z) cm will be at tem¬ 
peratures of hundreds to thousands of kelvin, approxi¬ 
mately 4 orders of magnitude greater than the 21 cm sig¬ 
nal. All the significant foreground contaminants should 
have smooth power-law spectra. Known sources of radio 
recombination lines are estimated t o contr i bute to the 
fluctuations at an insignificant level IQh fc MaddEnnft . 

Before fitting a model to the cosmological signal, it 
is necessary to clean the foregrounds from the data. 
The idea is to subtract out a smooth function from 
the total signal prio r to the parameter fitting stage 
llTegmark et al. 1 1200011 . Such pre-processin g is common 
with CMB data sets, and IWang et all 11200511 showed that 
this procedure can also be used in handling 21 cm obser¬ 
vations. 

At the frequency Vi in a pixel with angular index k±, 
an interferometer measures 

Xi = s(k±,Ui) + n(k ±1 Vi) + f(k±, z^), (26) 

where s is the 21 cm signal in visibilities, n is the detec¬ 
tor noise fluctuation, and / is the foreground amplitude 
(all of which are complex). We will subsequently write 

7 To more accurately capture these modes, we could discretize 
k and physically count the number of modes within the volume. 
However, our approximation is only inaccurate for small k. Fore¬ 
grounds will eliminate our ability to measure these long-wav eleng th 
modes such that a more precise treatment is unnecessary (. 

8 The reader may be familiar with expressions for the error that 
cont ain a factor of' yA/ N c . We do not have a factor of yA in equa¬ 
tion EH because each pixel has both a real an imaginary compo¬ 
nent. Since we only count pixels in the half space, this formulation 
is equivalent. 
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the quantities x, s, n and / measured at the N frequen¬ 
cies v = (zq, i/jv) with resolution Az/ as the vectors 
x,s,n , and /. There is one key difference between our 
calculation and that of lWang et al. . 1 200 51. Rather than 
subtracting from log(ir) a polynomial in log(z'), which is 
functionally very similar to the known foregrounds, we 
instead subtract a polynomial in v from x. While this dif¬ 
ference may require a higher order function to adequately 
fit the data, it also permits an analytic treatment. 

Fitting an order-n polynomial to the vector x is 
equivalent to projecting out the Legendre polynomi¬ 
als Po, Pi,..., P n , normalized such that Pi x (2v IB — 
1 )Pi 2 (2v/B - \)dv = S iu i 2 and P iti = Pi(2ui/B - 
1) \j Aza 9 Projecting out to order n, our cleaned signal 
is 

*= (i-E PlP ?) x = E p i p ? x (27) 

V 1=0 ) Z=n+1 


and a, n and / are defined in analogy to x. The covari¬ 
ance matrix for the cleaned signal is C = ( xx 4 ). Let 
us write II = jEn+i Pi Pf ■ We need to invert C to 
calculate SP'£ l T . Because C is singular, to invert C we 
use the trick C —> C + r] (E" =1 7*; Pf'j = C*, where 

r] is a large number. This met hod for inverting C does 
not lose information llTegmarklfl997D . In the basis of the 

Pi 

n 

C, = (ix^)+r 1 Y / Pi P J 

1 = 1 


= T J v 2 n + // t +5S t +r ? ^PzPf, (28) 

1 = 1 


where Tjv is the AT/v defined in equation cej except 
with the replacement to —» tk ± ((n rcJ) is diagonal in the 
chosen basis). When the detector noise dominates over 
the signal, the inverse of C* is 


c; 


i 


j_ _ // f 

Tn 2 Tn 2 (T n 2 + f) ’ 


(29) 


where I is the identity matrix. Here, we have dropped 
terms proportional to 1/rj. If the foregrounds can be 
cleaned well below the signal, the Fisher matrix for the 
21 cm power spectrum is 


Fk,k' = 


d 2 (log£) 


dP^ T (k) dPf T {k’) 


= tr 


dC 


dC 


dPf T {k) * dP£ T (k') 


(30) 

(31) 


where C is defined in equation m and only visibilities 
with the same k± are used. 

We want to constrain the parameters Pat (fcj_, k\ \). We 
can write PAr(k±, k\\) in terms of the signal s via a 
Fourier transform: 


(s« f )fcx = E w P Ar( k ±, k )PkPi- ( 32 ) 

k 


9 The formalism discussed in this section should apply to any 
complete set of orthogonal functions and not just Legendre poly¬ 
nomials. 


Note that k in and in PAr(k±,k) denotes the LOS 
component of k rather than the norm of k. Here the 
Fourier vector /x fc oc exp [iy/Bkv], where y is the length 
of the box, and w ~ X 2 B 2 /(A e x 2 y) (see equation d. 
It follows from equation © that 

dc 

dPf T {k) = w n Pk Pk nT = w Mfc ftfc- (33) 


For the k and k' at which the foregrounds can be cleaned 
well below the signal, the Fisher matrix is 



Since pixels with different k± are independent, we can 
combine the error in P^ T (k) from all pixels with the 
same (k, 9) as in (0 Therefore, if cleaning is successful, 
the combined error from the pixels in an annulus indexed 
by (k, 9) is 

5Pf T {ki,6) » ; 1 ,/fp^ -1 ) . (35) 

sjN c {ki,9) V V >i,i 

This equation is a good approximation for k at which 

Ml (sS t )/x fc f l^k ■ If vze use the approximate 

orthogonality of the fi k (note that the LOS k are sampled 
2n/y apart), this condition reduces to 

P%r(k ± ,k) » M * Mfc 2 = Q k± (k, n). (36) 

w (pI Ait) 

The larger the bandwidth, the higher order polynomial 
it should take to fit the data. To optimize the foreground 
removal procedure, the minimum n should be chosen such 
that the condition given in nu is satisfied. The larger 
the value of n, the more power will be removed from the 
21 cm signal. 

The formalism in this section can be easily generalized 
to include the situation in which the foregrounds are re¬ 
moved over a larger bandwidth than the bandwidth over 
which the 21 cm signal is extracted. Increasing the band¬ 
width over which the foreground removal is performed 
will improve an interferometer’s sensitivity to the cos¬ 
mological signal 1 H5.2I1 . 


4.1. Foreground Model 

The three major foreground contaminants are ex- 
tragalactic point sources, Galactic bremsstrahlung and 
Galactic synchrotron. The Galacti c synchrotron com - 
prises about 70% of the foreground ( iShaver et all19991) . 
but the ex tragalactic point sourc es may be the hardest 
to remove (iDi Matteo et alj[2002l) . Here we are not con¬ 
cerned with the overall amplitude of these foregrounds, 
since an interferometer cannot measure the k± = 0 
mode. 

To model the angular power spectrum of the Galactic 
synchrotron, we employ the function 

l2ClM ( * ^>^>2), (37) 


27r 


T. sy » = Al yn ( --- ) 

lo w lo V150 MHz/ 


n — Ao 


i°g(i 


( 38 ) 





























g syn = 2. 55, Aa syn = 0 .1, 0 = 2.5 and A * y * = 25 K 
ll Shaver et al.l Il999t iTegmark et aL 1200 0; Iwang et alJ 
120051) . The latter two values are extrapolated from 
30 GHz CMB observations. 

For the extragalacti c point sources we employ the 
IDi Matteo et ah 


galacti c point sou rces we em ploy tne 
Il2002i) model. [Gnedin fc Shaved i 21 )l) 1 : 


points out that Di Matteo et al I l200i ~ make very pes¬ 
simistic parameter choices for this model. As a re¬ 
sult, this model probably overestimates the contribution 
from extragalactic point sources. The extragalactic point 
source contribution has two components, a Poisson com¬ 
ponent and a clustering component. Bright sources can 
be removed from the map prior to the foreground fit¬ 
ting stage. Once bright sources are cleaned, the Poisson 
component is 


r T cnt /• pjl j\r 

c r = l dTj dC T>—, (39) 

where T cut is the minimum brightness temperature of the 
sources that can be cleaned and d 2 N / dTd( is the number 
of sources per unit brightness temperature at 150 MHz 
per ( - the spectral index of a source - per steradian. 

To model the clustering term, we assume that the spec¬ 
tral indexes of sources are spatially uncorrelated and set 
the correlation function of the extragalactic sources to 
be w(9) = [0/(0*)] - ' 8 , such that 

Cf ust K v 2 ) <x {If- 2 T eg {v x ) T eg (v 2 ), (40) 
where j3 = 0.85, 0* = 4' (IDi Matteo et al.lll200 2) and 



dT„, 


d( 


d 2 N 

dT Vo d( 


T V0 (~) 

vo 



(41) 


We model the probability distribution of the spectral 
index ( as a spatially constant Gaussian wit h stan dard 
devia tion 5( — 0.3 and mean </ = 2.8 llTegmark et, alJ 
l2000l) . We assume 4 sources sr 1 mJy 1 at 880 mJy 
and a power-law scalin g in flux with exponent —1.75 
iffiilviatteo et al.112002 :1. Furthermore, we take T cut = 
7 x Ti ns t where the instrumental sensitivity limit is 


21 


A 2 21 


inst 


sys 


Aant A e \J 2 to B 


(42) 


The values of S cut = 6.9xlO 5 [T cut /(1K)[ [i//(150MHz)] 2 
mJy for MWA, LOFAR, and the Square Kilometer Array 
(SKA) after 1000 hr of observations at 150 MHz with 
B = 6 MHz are listed in Tabled 

The foreground power is dominated by the Galactic 
synchrotron at most scales. Because of this, this fore¬ 
ground is the most difficult to remove from the 21 cm 
map. At l > 5000, the extragalactic point sources fluc¬ 
tuations start to become important. In this analysis, 
we ignore the contribution owing to Galactic and extra¬ 
galactic bremsstrahlung emission. The Galactic emission 
is expected to account for roughly 1 % of the contam ina- 
tion at the relevant frequencies ll Shaver et al .11199^) and 
contributes a negligible amount of power at all scales. 
While there is large uncertainty in the extragalactic 
bremsstrahlung, its contribu tion will also be minor at 
the relevant scales llSantos et al.il2005S) . 

With this model for foregrounds, we can calculate the 
experimental sensitivities using the formulas in the first 
part of this section if we note that {f{v i)* f{v 2 ))k ± =i/x = 


X 2 /A e Ci{vi, u 2 ), where the prefactor of X 2 /A e comes 
from f dudv\A„(u, v)\ 2 (see eqn. EU- 

We have made several simplifying assumptions for the 
form of the foregrounds. For example, extragalactic 
point sources will not exactly have a Gaussian distri¬ 
bution of spectral indexes and the frequency dependence 
of the foregrounds may be a function of l. While we 
anticipate that our simplifications will have a negligible 
effect on the overall foreground cleaning, this is a ques¬ 
tion that is b eyo nd the scope of our analysis. (See, e.g. 
iSantos et al.l ( 20051’) for a treatment of more complicated 
foreground models.) 

5. SENSITIVITY OF UPCOMING INTERFEROMETERS 

The MWA, LOFAR and SKA instruments are in var¬ 
ious stages of design planning. 10 In our calculations, we 
try to be faithful to the tentative design specifications 
for each facility and to make reasonable assumptions re¬ 
garding features of each array that have not been pub¬ 
licly specified. Tabled lists most of the parameters for 
these arrays that we use for our sensitivity calculation. 
Unle ss otherwise stated, the parameters we adop t come 
from iBowman et al . (2005a!) for MWA. Ide Vod J20041) 
and w ww.lofar.org for LOFAR, and lOa.rilli Kr. Raw lings 
(200l for SKA. 


5.1. Interferometers 

LOFAR will have 77 large “stations,” each of which 
combines the signal from thousands of dipole antennae 
to form a beam of « 10 deg 2 . Each station is also able 
to simultaneously image N p regions in the sky. We set 
N p = 4 in our calculations, but this number may be 
higher. The signal from these stations is then correlated 
to produce an image. In contrast, MWA will have 500 
correlated 4 m x 4 nr antenna panels, each with 16 dipoles. 
This amounts to a total collecting of 7000 m 2 at z = 8, or 
15% of the collecting area in the core of LOFAR. While 
correlating such a large number of panels is computation¬ 
ally challenging, this design gives MWA a larger field of 
view (FOV) than LOFAR (800 deg 2 ), which is an advan¬ 
tage for a statistical survey. 

The properties of SKA have not yet been finalized, 
and it is quite possible that the EOR science driver for 
SKA may form a distinct array from the other, higher- 
frequency drivers. In addition, the successes of MWA 
and LOFAR will likely influence the final design of SKA. 
The collecting area for SKA is projected to be roughly 
100 times larger than that of MWA. There are currently 
several competing designs for SKA’s antennae. At one 
extreme, SKA will have roughly 5000 smaller antennae 
(like a much larger MWA). At the other extreme, it will 
have fewer than 100 large antennae, each of which can 
simultaneously image several regions of the sky. For our 
calculation, we use the former extreme case, which makes 
it easier to have shorter baselines and to smoothly sample 
points in the u-v plane,- both of which are important 
considerations for EOR interferometers. We assume that 
the collecting area for SKA scales as A 2 , like a simple 
dipole, and is equal to 6 x 10 5 m 2 within the inner 6 km of 
the array for A = 21 (1 + 8) cm. This scaling is somewhat 

10 PAST is furthest along in construction, but it is not included 
in our analysis because detailed specifications are not publically 
available. PAST’s collecting area is comparable to that of MWA. 
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unrealistic, and the scaling of the collecting area will also 
depend on the spacing of the individual dipoles because 
the antennae will inevitably shadow each other at the 
longer wavelengths. 11 

The exact antenna distribution has not been decided 
for any of these instruments. For all three interferom¬ 
eters, we assume that the distribution of baselines is a 
smooth function. 12 The distribution of baselines in an ar¬ 
ray can substantially impact the sensitivity to the EOR 
signal. For MWA, we cal culate the sensitivities for an r~ 2 
antenna density profile llBowman et alJl2005al) . Specifi¬ 
cally, this distribution has a core with a physical covering 
fraction close to unity out to 20 nr before an r~ 2 falloff 
and a sharp cutoff at 750 m. The baselines are not as 
concentrated for the other two arrays. LOFAR will have 
an inner core within 1 km that has 25% of its antennae 
and an outer core with radius equal to 6 km with another 
25% of its antennae. For SKA we take (20%, 30%, 5%) 
of the antennae within (1, 6, 12 km). For SKA(LOFAR) 
we ignore the antennae outside 12(6) km for our calcula¬ 
tions. For simplification, we also assume that the density 
of the antennae is constant within each outer annulus for 
LOFAR and for SKA. However, we choose the inner 1 km 
region of both arrays to have a similar r -2 distribution to 
MWA, except with a wider core prior to the r~ 2 falloff in 
the differential covering fraction. The lower limit on the 
baseline length is approximately 4 m for MWA and 100 m 
for LOFAR, and we set this to be 10 m for SKA, which 
is approximately the physical diameter of the antennae 
panels. 

For these three arrays, the system temperature is dom¬ 
inated by the sky temperature. In our calculations, we 
set T sys = T s k y = 250 K at z = 6 , T SVR = 440 K at z = 8 
and T sys = 1000 K at 2 = 12 iHnwman et a,Uf2005a ). 
and we set B = 6 MHz bandwidth, which translates to a 
conformal distance of 100 Mpc at z = 8. In Annendixf/Vl 
we discuss how the choice of bandwidth can affect obser¬ 
vations. For the sensitivity calculations in this section, 
we chose observations that minimized the thermal noise, 

- which is the dominant source of noise on most scales, 

- by restricting the observation for each array to a single 
FOV. Finally, we set Au = 0.01 MHz for all of the ar¬ 
rays. While these arrays will have even better resolution 
than this, improved frequency resolution does not affect 
our results. 

5.2. Results 

Figures El El and ED show the pixel imaging capabil¬ 
ity and the statistical error in k 3 -Pat(4')/27t 2 , for MWA 
( dashed curves ), LOFAR (dash-dotted curves ), and SKA 
(solid cu rves) . For these figures, we use the parameters 
given in (EH and assume 1000 hr of observation over a 6 
MHz band and that the signal comes from the Universe 
when Xi = 0 and T s 3> Xcmb- For different ionization 
fractions, the signal can be both larger and smaller than 
the assumed signal. This is illustrated in FigureEJ where 
the thin solid curves represent the fiducial signal and the 
thin dashed curves represent the signal in the FZH04 
model for x t = 0.2, 0.55, and 0.75 for z = 12,8, and 6, 
respectively. 

11 This assumes that the low frequency part of SKA consists of 
dipoles. 

12 For LOFAR, which has far fewer antennae units than the 
other arrays, this assumption of continuity is fairly crude. 


Figure 21 plots the cumulative number of Fourier pix¬ 
els for wavenumbers less than k that have ratios of the 
RMS signal to the RMS detector noise that are greater 
than unity. In this plot, we do not include k < 2ir/y in 
the summation because, as we will show, the foreground 
removal procedure makes it unlikely that we can detect 
the cosmological signal for values of k smaller than the 
depth of the survey. Because MWA has a large FOV and 
is able to measure shorter baselines than the other inter¬ 
ferometers, it “images” a number of Fourier pixels com¬ 
parable to the number from LOFAR, despite having less 
collecting area. The sensitivity of these interferometers 
inevitably declines with redshift, with the detector noise 
in a pixel scaling roughly as T s k y y / A 2 /A e ~ (1 +z) 2 6 , as¬ 
suming that A e ex A 2 and that T s ^ y oc A 2 6 . Both LOFAR 
and MWA will have fewer than 1000 high signal-to-noise 
ratio (S/N) pixels at redshifts 8 and almost no pixels at 
higher redshifts. An SKA-class experiment will be re¬ 
quired to image modes with k > 0.1 Mpc -1 or z > 10. 

Observations of high redshift 21 cm emission are 
promising for cosmology in part because of the much 
larger number of Fourier modes that these observations 
can probe compared to other cosmological probes. These 
experiments can potentially probe scales larger than the 
Jeans length at all times during which the universe is 
neutral. The CMB, on the other hand, can only probe 
primordial fluctuations up to the Silk damping scale 
(/silk ~ 4000) from a single angular power spectrum. 
Currently, CMB experiments can image /^ ax < / 2 ilk in¬ 
dependent modes. The Wilkinson Microwave Anisotropy 
Probe (WMAP) is cosmic variance-limited for modes 
smaller than Z max ss 400, and Planck is cosmic variance- 
limited by 4nax ~ 3000. The number of modes that can 
be imaged by SKA in a 1000 hr observation at z = 6 in a 
6 MHz band is larger than the number of imaged modes 
for WMAP and significantly less than this number for 
Planck. A longer observation or a larger bandwidth will 
increase the number of modes that these interferometers 
can observe. 

The reason the S/N is generally smaller for 21 cm mea¬ 
surements than it is for measurements of the CMB is in 
part due to the bandwidth of these observations. Both 
cosmological probes are looking for fluctuations that are 
of order 10 -5 times that of the sky temperature. Un¬ 
fortunately, the number of independent samples of the 
sky temperature is proportional to the bandwidth, and 
CMB experiments have ~ 1000 times larger bandwidth 
since they observe at v ss 100 GHz. Therefore, CMB ex¬ 
periments can beat down their uncertainty in T s k y by an 
additional factor of (1000) 1 / 2 for an observation of the 
same duration. 

Figure 0 plots the fraction of pixels for a given value of 
k with a ratio of RMS signal to RMS noise greater than 
unity for MWA, LOFAR, and SKA. The vertical hatched 
line indicates 2-7r/y, and it is likely that foregrounds can 
be cleaned well enough only at scales rightward of this 
line. A larger fraction of LOFAR’s pixels than MWA’s 
pixels will have high S/N. Because MWA has a larger 
FOV, it still can detect a comparable number of high S/N 
pixels (see Fig. OJ. 13 SKA will have high S/N detections 

13 LOFAR should, however, have more success imaging a single 
quasar on the sky than MWA because it can image a larger fraction 
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TABLE 1 

The parameters that we adopt for MWA, LOFAR, and SKA. a 


Array 

Want 

A”ant A e (m 2 ) 

* = 6/8/12 

FOV 

(deg 2 at z=8) 

5 cut b 

(/Hy) 

min. base¬ 
line (m) 

Cost 
(10 6 $) 

MWA 

500 

4500/7000/9000 

7r 16 2 

180 

4 

~ 10 

LOFAR 

64 

(3.5/4.2/7.2) x 10 4 

4 x tt2.0 2 

30 

100 

~ 100 

SKA 

5000 

(3.6/6.0/12.5) x 10 s 

tt5.6 2 

2 

10 

~ 1000 


a These are the parameters for a central region of LOFAR and SKA and not the full 
array. We optimize the design for SKA for observations of the EOR, while keeping the 
current gross specifications for this array. 


b Values for 10 3 hr of observation with B = 6 MHz at 150 MHz. 


in almost all of its pixels up to k ~ 0.3 Mpc -1 . Figure 
0 illustrates that if foregrounds contaminate more large 
wavelength modes than is assumed, MWA and LOFAR 
can have substantially fewer high S/N pixels. Alterna¬ 
tively, a larger bandwidth will result in more high S/N 
pixels. 

Figure [S] compares the interferometers’ ability to sta¬ 
tistically constrain PAr{k), ignoring the effect that fore¬ 
grounds have on the sensitivity. Even though Pat(^) is 
not spherically symmetric, we spherically average Pat, 
as well as the errors, for the purpose of this plot. Because 
of this averaging, these interferometers will be slightly 
more sensitive to some modes than this plot implies. At z 
= 6, the trend is as expected: SKA is more sensitive than 
LOFAR and LOFAR is more sensitive than MWA. Still, 
LOFAR’s gains over MWA are not proportional to the 
square of the collecting area, as we might naively expect. 
At higher redshifts, LOFAR and MWA are comparably 
sensitive on most scales. We also plot the sensitivity of 
MWA at z = 8 for a flat distribution of antennae rather 
than the fiducial r~ 2 distribution of antennae. In this 
case, MWA is substantially less sensitive at all scales. 
This contrasts with angular power spectrum measure¬ 
ments, where a flat distribution of antennae is always 
more sensitive at larger k than a tapered distribution. 

If all the arrays had the same normalized distribu¬ 
tion of baselines, and if the error on the measurement 
of Pat(^) hr a Fourier pixel scales inversely with the 
square of the differential covering fraction for that k, as 
it does for the angular power spectrum, LOFAR should 
be many times more sensitive than MWA to a given pixel, 
- at least in the case where detector noise is the domi¬ 
nant source of noise. Because MWA observes many more 
independent Fourier cells owing to a larger survey vol¬ 
ume, for statistical detections, MWA should fare better 
even with the same distribution of baselines. However, 
the normalized distribution of baselines is not the same 
for these arrays. All these interferometers have a similar 
covering fraction in the very center, since the maximum 
covering fraction is unity. Since MWA stacks all its an¬ 
tennae in the core, it does not have any baselines that 
probe k± > 0.5 Mpc -1 at z = 8. Unlike MWA, only 
a quarter of LOFAR’s antennae are in its core region, 
leading to a smaller fraction of its baselines that can ob¬ 
serve modes with k± < 0.5 Mpc -1 . Another reason that 

of its pixels. Since MWA is extremely cored, this will make its beam 
much more coarse than the other experiments. 


LOFAR is not many times more sensitive than MWA is 
because LOFAR’s minimum baseline of 100 m does not 
allow it to detect modes with k± < 0.03 Mpc -1 (Fig. |SJ). 
These modes happen to be those to which MWA is most 
sensitive. 

Figures [3 and 0 illustrate how foregrounds affect the 
sensitivity of MWA to the power spectrum for a fixed 
k i and to statistical detections of the power spectrum. 
The foregrounds are first fitted in the frequency direc¬ 
tion for each k±. As we will see, the foreground pre¬ 
processing removes significant power from the signal on 
large scales, reducing the sensitivity to the signal at such 
scales. In Figure [3 we assume a l-er angular fluctuation 
for the foreground model outlined in We then sub¬ 
tract a quadratic (left panel ) or cubic polynomial (right 
panel) from the foregrounds over a frequency interval of 
6 , 12, and 24 MHz. For these figures, the 6 MHz band 
from which we extract the 21 cm signal is centered in the 
larger frequency intervals in which we remove the fore¬ 
grounds. The placement of this band does not affect the 
results substantially. We find that for all the bandwidths, 
a quadratic or a cubic polynomial is able to remove the 
residual foregrounds Qk ± , defined in equation sub¬ 
stantially below the signal (thin solid line). The solid, 
dashed and dash-dotted curves at the bottom of Figure 
0 indicate Qk± for foreground removal in a 6, 12 and 
24 MHz band. The cubic polynomial is able to remove 
the foregrounds well below the signal for all cases, but 
for B < 12 MHz (or at larger k± than shown, where the 
foreground contamination is smaller) the quadratic poly¬ 
nomial is sufficient. This conclusion holds for the other 
interferometers as well. 

Foreground subtraction removes power from the cos¬ 
mological signal on smaller scales as we decrease the 
bandwidth over which we remove the foreground (or as 
we increase the order of the polynomial). Our analysis 
accounts for this by effectively dividing the sensitivity 
curves by the filter function that describes how power 
is removed from the cosmological signal as a function 
of fc|| (Fig. 0 upper thick curves). These sensitivity 
curves would be flat in the absence of foregrounds and 
foreground cleaning. The foreground cleaning causes 
these curves to be less sensitive to the signal at large 
scales. The errors on the power spectrum are substan¬ 
tially smaller for the second generation of interferome¬ 
ters, but the foreground removal has a qualitatively sim¬ 
ilar effect on these interferometers’ sensitivity curves. 
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We can also combine the signal along different k± in 
Fourier space to constrain P\r(fc) (eqn. EEJl- Figure 
0 shows the statistical error in the spherically averaged 
21 cm power spectrum for MWA after foregrounds are 
cleaned. The thick solid line represents the error just 
from detector noise. At around the scale correspond¬ 
ing to the depth of the box, too much power from the 
cosmological signal is removed due to foreground clean¬ 
ing for MWA to be sensitive. As before, this effect is 
minimized by fitting to a larger bandwidth. The dashed 
and dash-dotted curves indicate the errors if we remove 
the foregrounds with a cubic polynomial (n = 3) in a 6 
and 24 MHz band, respectively. The thick dashed curve 
shows the errors for a quadratic polynomial (n = 2 ) with 
B = 6 MHz: n = 2 removes substantially less power 
than n = 3 given the same bandwidth. Similar conclu¬ 
sions hold for other EOR interferometers. 

Despite the simple model we employ for foregrounds, 
we expect our conclusions pertaining to foreground re¬ 
moval to be fairly robust. Our technique should be able 
to remove the foregrounds from a mode with fc|| 2ir/y 
as long as P? (k) <C P 21 (fc), where Pf is the foreground 
3-D power spectrum. The foregrounds are expected to 
be in this limit for most relevant k. The process of re¬ 
moving the foregrounds from the signal will inevitably 
remove the signal for k\\ < 2ir/y. 

This is not to say that removing foregrounds from 21 
cm maps is trivial. Our analysis neglected several compli¬ 
cations that the real observations must deal with. Since 
the observed wavelength increases with redshift, over the 
depth of the survey a mode with a set value of k± will be 
measured by different baselines. In our analysis, we ig¬ 
nored this effect. As long as the distribution of baselines 
is fairly smooth, we expect that this will have a minor 
effect on foreground removal. Other foregrounds that are 
beyond the scope of this paper include residuals owing 
to imperfect point source subtraction, radio frequency in¬ 
terference contaminating frequency intervals within the 
observation band (this may be a substantial challenge for 
LOFAR, which is in a radio loud environment), and the 
residuals that arise owing to the imperfect modeling of 
atmospheric distortions (modeling the atmosphere may 
be a significant challenge for MWA due to its large FOV). 

6. COSMOLOGY FROM THE 21 CM POWER SPECTRUM 

Observations of high redshift 21 cm emission are ca¬ 
pable of measuring Pss on smaller scales than current 
CMB experiments. Our calculations show that SKA can 
sensitively probe comoving megaparsec scales, which are 
also smaller than scales observed by galaxy surveys and 
comparable to scales probed with the Lya forest. The 
sensitivity to smaller scales than the CMB may allow 21 
cm observations to break degeneracies among cosmolog¬ 
ical parameters that are present in CMB constraints. 

In this section, we utilize the sensitivity calculation 
described in m and $2 to estimate how well upcoming 
EOR interferometers can constrain cosmological param¬ 
eters from Pat, and, in particular, whether these con¬ 
straints will be competitive with CMB observations. We 
divide our discussion into three cases: (1) If density fluc¬ 
tuations dominate the signal (ED- This can happen if 
Xi <C 1 and Tomb T s or if X-rays are responsible for 
the reionization of the Universe. (2) When the bubbles 
contaminate the P^o and P ^2 terms such that only P M 4 



Fig. 4. — Integrated number of Fourier pixels Ngn v>i(^) with 
k' < k that have a ratio of the RMS signal to the RMS detector 
noise that is greater than unity for a 1000 hr observation with 
B = 6 MHz. We use the specifications given in Tabic HI and in 
33 for MWA ( dashed curve), LOFAR (dot-dashed curve), and SKA 
(solid curve). These curves do not include pixels with k < 2n/y, 
since foregrounds will contaminate these pixels substantially. The 
21 cm signal for this calculation is from a fully neutral medium in 
which T s Tq mb- If the universe is partially ionized at 2 : = 8, 
the signal can be both larger and smaller than this (see Fig. |6j. 
By a redshift of 12, only SKA will have any high S/N pixels. 



Fig. 5. — Fraction of Fourier pixels for MWA ( dashed curve), 
LOFAR (dot-dashed curve), and SKA (solid curve) that are 
“imaged”— that is, they have a ratio of RMS signal (assuming 
a neutral universe) to RMS detector noise that is greater than 
unity—after 1000 hr of observation in a Fourier shell of radius k. 
The hatched vertical line marks the depth of this 6 MHz observa¬ 
tion at z = 8. Scales to the left of this should be wiped out by 
foregrounds. LOFAR can image a substantially higher fraction of 
pixels at the relevant k than MWA, and SKA can image nearly all 
of its pixels up to k = 0.3Mpc —1 . 


and - which arises from the AP effect, - are pris¬ 
tine enough to measure cosmological parameters (E3- 
(3) On large scales at which neutral fraction fluctuations 
are important, but at which these fluctuations trace the 
density fluctuations (H3J). Note that the analysis in this 
section does not assume any mo del for reionizat i on. 

It came to our attention that (Bowman et al.1 (l2005bfl 



















Fig. 6. — Detector noise plus sample variance errors for a 1000 hr 
observation on a single field in the sky, assuming perfect foreground 
removal, for MWA (thick dashed curve), LOFAR (thick dot-dashed 
curve), and SKA (thick solid curve), using the specifications given 
in Table and in 33 and for bin sizes of A k = 0.5 A:. These errors 
are for the spherically averaged signal (see text for the discussion of 
this point). The hatched line in the middle panel represents MWA 
with a flat distribution of antennae rather than the fiducial r ~ 2 
distribution. The detector noise dominates over sample variance 
for these sensitivity curves on almost all scales. The thin solid 
curve represents the spherically averaged signal for Xi 1 and 
T s Tqmb- We use this curve to calculate the sample variance 
error. For comparison, the thin dashed curves show the signal 
from the FZH04 model when Xi is equal to 0.20, 0.55 and 0.75 for 
z = 12,8 and 6, respectively. 


was performing a similar analysis for MWA and 
MWA5000. This paper, submitted concurrently with 
ours, pertains to when the signal is in the regime we 
discuss in a There are a few differences between our 
two approaches . In add ition to cosmological parameters, 
iBowman et al ] (I2005H1 fits to parameters for the fore¬ 
groun d residuals as well as o ther observational parame¬ 
ters. IBowman et alJ ( 2005b ) also assumes a spherically 
symmetric Pat- Deviations from spherical symmetry en¬ 
hance cosmological parameter constraints. Another im¬ 
portant difference is that our analysis combines 21 cm 
observations with current and future CMB experiments. 
This can break parameter degeneracies present in these 
separate cosmological probes, and is important for assess¬ 
ing the true value of 21 cm observations for cosmology. 

6 .1. Density Fluctuations Dominate 

We first concentrate on the signal in the case where 
the density fluctuations dominate over the spin tem¬ 
perature and neutral fraction fluctuations. This is the 
case in which 21 cm observations will be most sensi¬ 
tive to cosmological parameters. If reionization occ urred 
at z sa 6.5, as the Sl oan quasars suggest lIBecker et all 
12001 iFan et al.H20f)l . it is not altogether unlikely that 
upcoming interferometers will observe the signal in this 
regime. Models in which x-i = Q f co ii have a significant 
period in which ij < 1. During this period, fluctua¬ 
tions in xh are unimportant. In addition, it is expected 
that at higher redshifts than those considered here, X- 


Quadratic Cubic 



Fig. 7. — Foreground removal for MWA for a 1000 hr observation 
with B = 6 MHz. The thin solid curve shows Pat{^±, h\\ ) plotted 
as a function of fc| | for a neutral universe. The thick upper curves 
show the sensitivity to the signal <5 Pat (fc \ -T| |) if we remove the 
foregrounds using the method outlined in 31 in a band centered 
on the 6 MHz cosmological signal processing window of 6, 12, or 
24 MHz (solid, dashed or dot-dashed curves, respectively). Note 
that the sensitivity curves are above the curves for signal because 
this is the sensitivity to a single Fourier pixel ( N c = 1). The 
lower curves show the function Qk ± , ~~ the residual foreground 
level, - for each of the three processing window bandwidths (see 
33 ). The foreground power in each panel is assumed to be a 1-cr 
fluctuation of the model discussed in 31 The left panel shows a 
fit with a quadratic polynomial (n = 2) and the right shows a fit 
with a cubic polynomial (n = 3). While increasing n will always 
remove more of the foregrounds: it will also remove more of the 
cosmological signal as well. This can be seen in this plot by the 
reduced sensitivity to the signal at large scales. (In the absence of 
foregrounds and foreground cleaning, the sensitivity curves would 
be flat in this plot.) 


ray photons and shocks heat th e gas in the IGM to well 
above the CMB temperature ( Venkatesar^^n J2CKBJ; 
IChen fc Mi r al d a- Es c ud el 12 0 0 -111 . Cjanj^^Maxjanl II2003I 
argue that around z = 20 the first stars will produce 
a large enough background in Lyct such that T s will be 
coupled to the kinetic temperature of the gas through 
the Wouthuysen-Field effect. If this is true, spin temper¬ 
ature fluctuations will be subdominant at the redshifts 
we consider. 

Tables |21 and |o] quantify how the 21 cm signal can con¬ 
strain some of the most interesting cosmological param¬ 
eters: r, 0, w ,w, D, m h 2 , Di,h 2 , n s , Sjj, ct s and The 
tilt n s we define to be the power law index of the pri¬ 
mordial power spectrum at k = 0.05 Mpc^ 1 and a s = 
dn s {k)/d\og{k). The parameter Sh is roughly the size of 
density fluctuations at the present day horizon scale (de¬ 
fined here such that the primordial power spectrum to¬ 
day P(k) = 27r 2 <5^fc n ‘»( fe )/[70kms- 1 Mpc- 1 ] 3+n »( fe )). To 
construct the linear power spectrum used in this analysis, 
we employ the transfer function from the code CAMB. 14 
To get c onfidence intervals, w e use the Fisher matrix for¬ 
malism ( Tegmark et a l. (l997f) . 

In the long term, MWA plans to increase the num¬ 
ber of antenna panels from 500 to 5000. This array, 

14 Available at http://camb.info 
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Fig. 8. — Dimensionless 21 cm power spectrum for a neutral 
universe measured in a 6 MHz band (thin solid curves) and the 
1-cr error for a 1000 hr observation by MWA in the absence of 
foregrounds (thick solid curves ), after foreground subtraction with 
a cubic polynomial at all k± in a 6 MHz window (dashed curve) 
and in a 24 MHz window (dot-dashed curve). The thick dashed 
curves is for a quadratic polynomial from a 6 MHz band. The ver¬ 
tical hatched line denotes the scale of the 6 MHz box. Foreground 
cleaning reduces the sensitivity at large scales. This plot illustrates 
that it makes a big difference how foregrounds are removed from 
the signal: experiments will want to remove the foregrounds over 
a fairly large bandwidth and with as low order a polynomial as 
possible. 


MWA5000, will have a comparable collecting area to LO- 
FAR and 10% of the collecting area of SKA. To model 
MWA5000, we use an r~ 2 distribution of antennae out to 

1 km, similar to MWA, but with a larger flat core than 
MWA that extends out to 80 m rather than 20 m. We 
also include MWA50K, which is another 10 times larger 
than MWA5000, but again built in the same mold as 
MWA. Correlating 50,000 antennae will be a significant 
computational challenge. 

In Table |21 we calculate the 1 -a errors on cosmological 
parameters for observations at z = 8 in the case in which 
Xi < Cl and T s Xcmb, such that all the /r terms trace 
the density power spectrum. Unless otherwise noted, we 
perform the calculations in this section for an observa¬ 
tion of 2000 hr on two locations in the sky, or roughly 

2 productive years. 15 Observations should generally be 
chosen to minimize the number of patches on the sky 
and maximize the duration for which each patch is ob¬ 
served because detector noise dominates the uncertainty 
at most scales. For the second generation of EOR in¬ 
terferometers, this is not necessarily the case, and some¬ 
times parameter estimates are improved by choosing a 
different observing strategy. 

Future observations have the potential to improve 
many of the current constraints on cosmological param¬ 
eters. Two years of observation with MWA and LOFAR 
have trouble constraining a five parameter cosmology: 
Ua, flbh 2 , n s , and the normalization parameter 

xh Sh (see Table 0 note ‘b’). However, when combined 

15 Interferometric observations have never been integrated for 
such long periods on a single field. It is uncertain whether such 
observations are even possible, and this will depend heavily on how 
well we can deal with various systematics in these systems. 


with current CMB observations (WMAP, Boomerang, 
ACBAR [Arcminute Cosmology Bolometer Array Re¬ 
ceiver], and CBI [Cosmic Background Imager]), both 
measurements by MWA and LOFAR are able to improve 
measurements of Cl m h 2 , f2„, n Sl and a s . MWA and 
LOFAR are not able to significantly improve the con¬ 
straints from Planck. Unfortunately, Xh is not well con¬ 
strained by measurements by MWA or by LOFAR when 
they are combined with current CMB observations. This 
is because the current uncertainty in cosmological param¬ 
eters leads to substantial uncertainty in the amplitude of 
Ps L s L at relevant scales. Planck will be able to refine the 
measurement of these parameters and the first genera¬ 
tion of 21 cm experiments plus Planck will place tighter 
constraints on xh- 

The second or third generation of 21 cm observations 
will be substantially more sensitive to the cosmology. 
By themselves, MWA5000, SKA and MWA50K can con¬ 
strain a seven parameter cosmology that involves a s and 
Q. v in addition to the other five parameters we used for 
LOFAR and MWA (Tabled• Surprisingly, MWA5000 is 
comparably sensitive to SKA despite having 10 times less 
collecting area. This is because SKA is not as centrally 
concentrated with only 20% of its antennae in the 1 km 
core while MWA5000 has 100%. Also, MWA5000 has 
a larger FOV than SKA, which results in smaller errors 
on large scales, scales at which these arrays are sample 
variance limited. Large scales probe the baryonic wiggles 
and therefore can provide substantial constraining power 
(Fig. 0 . If we alter the observation for SKA to decrease 
the sample variance, - having it observe ten locations 
on the sky each for 400 hr, - then SKA’s sensitivity is 
improved (see entries with asterisks in Table 0 . 

In combination with Planck, MWA5000 and SKA can 
improve constraints on f l w , £l m h 2 , n s ,a s and fU, and 
MWA50K can do even better (Table 0 . Because these 
observations probe smaller scales than the CMB, the pa¬ 
rameters that affect the small scale behavior, namely 
n s ,a s and f2„, show the most substantial improvement. 
In addition, as one changes the cosmological parameters 
in the conversion from u to fc, this distorts the measured 
power spectrum in k-space, providing an additional ef¬ 
fect that can be used to constrain parameters. This is 
illustrated in Table 0 the constraints on SKA denoted 
with superscript ‘c’ are if we do not vary cosmological 
parameters in the conversion from u to k. The uncer¬ 
tainty on Ua, U m /i 2 , and £l v is substantially larger in 
this case. 

Tabic 0 shows the sensitivity to cosmological parame¬ 
ters if the above scenario occurs at higher redshifts than 
z = 8. MWA5000 is still sensitive to the signal at z = 10, 
and its sensitivity falls off at z = 12. Our design for SKA 
is unrealistically optimized for all considered redshifts. 
Because of this, SKA is more sensitive than MWA5000 
at z = 12, but it sensitivity is still falling due to the 
increasing sky temperature. 

The uncertainty estimates in this section are for obser¬ 
vations with B = 6 MHz. Experiments will be able to 
process a much larger bandwidth, and, if we are fortu¬ 
nate, nature could provide an even larger redshift slice 
in which density fluctuations dominate. 

6 .2. Neutral Fraction Fluctuations are Significant 
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TABLE 2 

Errors on cosmological parameter estimates when density fluctuations dominate the 21cm signal. 3. 



T 

Ok; 

W 


h? 

n s 

5 h x 10 s b 

a s 


X H 


0.1 

0.7 

-1.0 

0.14 

0.022 

1.0 

3.91 

0.0 

0.0 

1.0 

LOFAR 

- 

0.07 

- 

0.11 

0.03 

0.11 

5.0 

- 

- 

- 

MWA 

- 

0.06 

- 

0.09 

0.02 

0.09 

4.2 

- 

- 

- 

MWA5000 

- 

0.005 

- 

0.008 

0.002 

0.03 

0.37 

0.010 

0.007 

- 

SKA 

- 

0.005 

- 

0.009 

0.002 

0.06 

0.51 

0.016 

0.015 

- 

SKA C 

- 

0.11 

- 

0.042 

0.003 

0.07 

2.0 

0.017 

0.08 

- 

SKA* * 

- 

0.004 

- 

0.007 

0.002 

0.03 

0.32 

0.010 

0.008 

- 

MWA50K* 

- 

0.002 

- 

0.004 

0.001 

0.01 

0.17 

0.004 

0.002 

- 

CCMB 

0.060 

0.084 

- 

0.017 

0.0014 

0.072 

0.29 

0.039 

0.12 

- 

CCMB+ LOFAR 

0.057 

0.050 

- 

0.010 

0.0012 

0.027 

0.22 

0.022 

0.02 

0.2 

CCMB+ MWA 

0.056 

0.046 

- 

0.009 

0.0011 

0.021 

0.22 

0.022 

0.02 

0.2 

CCMB+ MWA5000 

0.048 

0.005 

- 

0.003 

0.0009 

0.013 

0.18 

0.005 

0.004 

0.06 

CCMB + SKA 

0.048 

0.005 

- 

0.003 

0.0009 

0.014 

0.18 

0.005 

0.007 

0.06 

Planck 

0.0050 

0.029 

0.09 

0.0023 

0.00018 

0.0047 

0.026 

0.008 

0.010 

- 

Planck +MWA5000 

0.0046 

0.017 

0.06 

0.0009 

0.00012 

0.0033 

0.018 

0.003 

0.003 

0.03 

Planck + SKA 

0.0046 

0.021 

0.08 

0.0008 

0.00012 

0.0034 

0.018 

0.003 

0.004 

0.04 

Planck + SKA* 

0.0046 

0.017 

0.07 

0.0007 

0.00012 

0.0032 

0.017 

0.003 

0.003 

0.03 

Planck + MWA50K* 

0.0045 

0.007 

0.03 

0.0004 

0.00010 

0.0029 

0.016 

0.002 

0.001 

0.01 


a We assume, unless otherwise noted, observations of 2000 hr on two places in the sky in a 6 MHz band which 
is centered at z = 8. Current CMB (CCMB) are the combined results of WMAP, Boomerang, ACBAR, and 
CBI. In these calculations, we account for foregrounds by imposing a sharp cutoff in sensitivity at k = 2n/y, 
where y is the width of the box, and we avoid fitting to scales in the non-linear regime by imposing a small scale 
cutoff at k = 2 Mpc - . These calculations are for a flat universe, 1 = and dashes indicate parameters 

which are not marginalized. 

b From just the 21 cm data, the parameter Sh is completely degenerate with xh • Because of this, for 21 cm 
observations alone, the constraints in this column are really for the parameter xh Sh- 

c We use the fiducial cosmology in the conversion from u to k such that the angular diameter distance and the 
depth of the map do not change when we vary parameters to get the above confidence intervals. 


* Observations of 10 locations on the sky, at 400 hr each. 


TABLE 3 

Same AS TABLED BUT FOR HIGHER REDSHIFTS. a 



T 

S2a 


Qf, h 2 

n s 

5 h x 10 5 b 

&S 


X H 


0.1 

0.7 

0.14 

0.022 

1.0 

3.91 

0.0 

0.0 

1.0 

MWA5000 (z = 10) 

- 

0.010 

0.014 

0.004 

0.03 

0.6 

0.01 

0.010 

- 

SKA (z = 10)* 

- 

0.007 

0.010 

0.003 

0.03 

0.4 

0.01 

0.009 

- 

MWA5000 (z = 12) 

- 

0.019 

0.030 

0.008 

0.07 

1.4 

0.03 

0.016 

- 

SKA (z = 12)* 

- 

0.011 

0.014 

0.004 

0.05 

0.7 

0.02 

0.013 

- 

Planck 

0.0049 

0.011 

0.0023 

0.00017 

0.0047 

0.03 

0.007 

0.010 

- 

Planck + MWA5000 (z = 10) 

0.0047 

0.007 

0.0013 

0.00013 

0.0036 

0.02 

0.005 

0.003 

0.03 

Planck + SKA (z = 10)* 

0.0046 

0.006 

0.0011 

0.00013 

0.0035 

0.02 

0.004 

0.003 

0.03 

Planck + MWA5000 (z = 12) 

0.0049 

0.009 

0.0017 

0.00015 

0.0040 

0.02 

0.007 

0.004 

0.04 

Planck + SKA (z = 12)* 

0.0047 

0.007 

0.0014 

0.00014 

0.0037 

0.02 

0.005 

0.004 

0.04 


a See table note (a) in Table [21 for the specifications used in these calculations. 

b From just the 21 cm data, the parameter Sh is completely degenerate with xh- Because of this, for the 21 cm 
observations alone, the constraints in this column are really for the parameter xh Sh ■ 

* Observations of 10 locations on the sky, at 400 hr each. 


In this section, we investigate whether it is possible 
to extract P M 4 = x 2 H Ps L s L from PAr(fc) well enough to 
constrain cosmological parameters when ionized fraction 
fluctuations are important. 16 Figure ITT)! shows the 2 : = 

16 The techniques in this section also apply to periods during 
which spin temperature fluctuations are important. 


8 sensitivity curves for MWA and for SKA. The right 
panel is shown at the beginning of reionization (x* = 0 . 1 ), 
when the density fluctuations are still the largest source 
of fluctuations. In this case, SKA will be sensitive to P M 4 
over 1-2 decades in k. Conversely, MWA is not sensitive 
to P^4. The left panel shows the opposite case, when 
the bubbles dominate (x, = 0.7). In this case, MWA and 
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k (Mpc -1 ) 


Fig. 9. — Cosmic variance plus detector noise errors on the 
large scale modes of the spherically averaged power spectrum for 
SKA ( thick error bars) and MWA5000 ( thin error bars) after 1000 
hr of observation in a 6 MHz band. We accentuate the wiggles 
in Pat W by subtracting a power spectrum that does not have 
baryonic wiggles P^,(k). The vertical hatched line indicates the 
size of a 6 MHz box, approximately where we expect foregrounds to 
swamp the signal. The second generation of interferometers will be 
able to detect the wiggles at a few a significance level. MWA5000 
is more sensitive than SKA to these features because of its larger 

POV TVio dicrnnl Vi<=>r£> i<a fnr 'v ■ 1 • +.Vi£» r>f TTTT VyiiVJiIoq 

□ 



Fig. 10. — Signal using the FZH04 model ( thin curves) and 
detector noise plus cosmic variance errors ( all other curves) at z = 
8 for a 2000 hr observation combining two different fields of view, 
or roughly 2 yr of observations, for the P o ( solid curves) and 
4 (dashed curves) components of the signal. In the left panel 
xi = 0.7, and in the right panel Xi = 0.1. We plot the sensitivity 
curves for MWA (medium-width curves) and for SKA ( thick curves) 
with a binning width of A k = 0.5 k. These curves assume that the 
angular diameter distance to z and that the value of H (z) are well 
constrained such that we can convert between u and k. If this is 
not the case, there is also a /A term. 


SKA are both not sensitive to P^4 . Both MWA and SKA 
can be sensitive to the P^o over a range of scales. This 
analysis assumes that the fiducial cosmology is correct or 
else we could not measure P^o and P^4, since we need to 
know the angular diameter distance and H(z) to be able 
to convert u to k. 


If the cosmology assumed in the conversion from u to 
k is incorrect, there will be an asymmetry in the rnea- 
su red depth versus measu red angular size of an object 
llAlcock fc Paczvnskl 11979 1. As a result, features, such 
as the effective bubble size, will appear distorted in the 
21 cm map [the AP effect]. Tests for this asymmetry 
can be used to put constraints on cosmological parame¬ 
ters. This effect is implicitly included in the analysis in 
flfi.1l because we vary the cosmological parameters in the 
conversion from u to k. To extract cosmology from P M 4 
alone, it is more subtle to include this effect. 

The presence of this new angular dependence from the 
AP effect will contribute an additional term to the p de¬ 
composition of Pat (equation 0, which has a unique 
U 6 angu lar dependence with coefficient |n usse 3 [20051 


Pu?{k) = -ax H ( 4 P SlSl - 


dP Sl 


Sl 


(9 log k 


), 


(43) 


keeping terms to linear order in a. Here 1 + a is the ratio 
of the true value of the Hubble constant times the angu¬ 
lar diameter distance, H(z) Da(z), to the value assumed 
in the conversion from it to fc. This term arises from 
the coupling of the AP effect to P^i. If detected, the p 6 
term would indicate a clear problem with the assumed 
cosmological model. The P^o, P ^2 and P^i terms are 
also affected by the AP effect. Since we currently can¬ 
not model P^o and P M 2 when the bubbles are important, 
the AP effect will not be detectable from these terms. 
However, we can model P^4 to zeroth order in a. If the 
assumed cosmol ogy is i ncorre ct such that a 7^ 0, this 
term becomes llBarkanall2005l l 


P M 4 = Pfi + a (5 P‘ r 4 - 2 P M 2 + 


/ dP u i 

-ol(3 P,5+ m 


d log k 


dP^ 

d log k 


(44) 


where P*4 is the true /i 4 term (what we measure if a = 0) 
and (1 + aj_) is the ratio of the assumed angular diam¬ 
eter distance Da{z ) to its true value. Since the P M 2 is 
generally larger than P*4, the deviation from P*4 can be 
quite significant. 

In Figure El we plot P* 4 ( solid curves) , P^e (dot- 
dashed curves ), and P^4 — P*J (dashed curves) for Xi = 
0.2 (thick curves) and Xi = 0.55 (thin curves) using the 
FZH04 analytic model at 2 = 8. We take a = 0.1 and 
ax = 0.1 to roughly match the current uncertainty in 
these parameters. The first generation of 21 cm arrays 
will not be sensitive to the AP effect or to P /J 4. The 
dashed and dot-dashed curves in Figure HT1 that are la¬ 
beled “Errors” represent the sensitivity curves for SKA 
to the P M 4 and P^ terms assuming 2 yr of observation. 
SKA will not have the sensitivity to give a meaningful 
measurement of either a or P M 4, but, after a sufficiently 
long integration, may be able to detect P^e. A similar 
conclusion holds for MWA5000. 

Let us now make this analysis more quantitative. If we 
take the Fisher matrix Fij at a given k for the parame¬ 
ters P M o, P^2, P^4, and P M s (indexed 1 through 4), then 
marginalize the contaminated parameters P^o and P M 2, 
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the Fisher matrix of just the parameters P^i and P^e is 



By the chain rule, it follows that the Fisher matrix for 
the cosmological parameters Ai,A„ obtained from just 
4 and P^e (indexed by 1 and 2 ) is 


F ij( k ) ~ F n ~g\~~d\ ; 

dP^ OP 

+ *21- 


dXi d\ 


dP^ 4 dP^e 

12 dXi dXj 
dP > dP 
22 dXi dX 3 ' 


(46) 


In Table 01 we consider the scenario in which xh = 0.8 
and only P^i and P^e yield a pristine measure of the 
linear density power spectrum via equations (TH 1 ) and 
(03) . (While P^i depends on P^ to linear order in a, 
since we can measure 2 , we can always use P ^4 and 
P M e to measure Pss-) We assume that a and aj_ are 
zero, but nonzero values do not change the results signif¬ 
icantly. In this case, MWA5000 and SKA can improve 
constraints moderately on cosmological parameters ob¬ 
tained from current CMB data sets. However, they are 
unable to compete with Planck. The measurement of 
P M 4 and P^e , when combined with Planck, will be most 
useful for measuring xh rather than for constraining cos¬ 
mological models. Two years of observation with SKA 
can constrain x h to better than 10%. Even MWA50K 
is not able to sensitively constrain the signal on its own 
with this method. 

SKA does fare better than MWA5000 for the analy¬ 
sis in this section, which was not the case in K6.1 1 This 
stems from MWA5000 being significantly more concen¬ 
trated than SKA and therefore not as sensitive to small 
scale modes perpendicular to the LOS. These modes 
are important to be able to separate the different p 
terms. Our design for MWA50K is very concentrated, 
like MWA5000, and so is also not optimal for measuring 
P M 4 and P^e. 

If one assumes some simple parameterization for P xx 
and P x s (or P M 0 and P M 2 ), rather than marginalizing over 
P M 0 and P ^2 for each fc-bin as is done here, then perhaps 
one could measure Pss with more confidence. 


6.3. Large Scales 

Up until now, we have ignored all components of the 
signal that are contaminated by the bubbles. On large 
scales, this may not be necessary. On scales much larger 
than the effective HII bubble size P e ff, if the Poisson 
fluctuations owing to the bubbles are unimportant, the 
bubble fluctuations will trace the density fluctuations. 
Thus, when k <C PA 1 , we may have the relations 

Pxx{k ) ~ b\ Pss, P x s(k ) « b 1 Pss■ (47) 

It is not necessarily the case with HII bubbles during 
reionization, as it is typically for galaxy surveys, that 
Poisson fluctuations are unimportant at large scales. If 
we include both the part of the signal that traces the 
Pss and the Poisson component, we can write the 21 cm 
power spectrum at large scales as 

Pat(^) = T 2 { \pH + ^1 ~~ ‘Zxh 61] + 2/i 2 

x [x 2 h - x H 61 ] + p 4 [x 2 h ] }P S s + P P oi,(48) 



Fig. 11.- Plot of k 3 P(k)/2n 2 for P(k ) equal to P *4 (solid 
curves), P^s (dot-dashed curves) and [(,4 — P l1 4 (dashed curves) 

for Xi = 0.2 (thick curves) and Xi = 0.55 (thin curves) at z = 8 . 
Here we assume that a = 0.1 and a± = 0.1. The dashed and 
dot-dashed curves, labeled “Errors,” are SKA’s errors for a 2 yr 
observation on P j: 4 and P^e when Xi = 0 . 2 , computed with bins 
A k = 0.5 k. 


where P po i is the Poisson contribution, which is constant 
in k. This type of decomposition may also hold when spin 
temperature fluctuations are important. On large scales 
we may again parameterize spin temperature fluctuations 
with the relation 5t s = &t s 5 for some constant bx 3 , and 
again the 21 cm power spectrum is proportional to Pss- 
If we include the AP effect, terms enter with derivatives 
of Pss- 

Equation 03 is promising in that, when Poisson fluc¬ 
tuations are unimportant, there are effectively three un¬ 
knowns (bi,x-H and Pss) and three equations, since the 
p 0 , p 2 and fi 4 components can, in principle, be measured. 
All three of the unknowns are very informative: Xh{z) 
tells us about the global reionization history, bi(z) indi¬ 
cates where the bubbles are located (i.e. overdense or 
average density regions), and of course Pss is quite in¬ 
teresting. In addition, b\ can be much larger than unity, 
such that the signal is enhanced over that of a neutral 
medium. 

There are several difficulties with extracting cosmol¬ 
ogy from these large scale modes. One difficulty is that 
modes larger than the width of the 21 cm survey will be 
contaminated by foregrounds. Another complication is 
that 21 cm interferometers will have progressively more 
trouble capturing terms of increasing order in n, which 
is necessary to separate the terms in equation 03 - De¬ 
spite these difficulties, it is probable that modes on scales 
near the baryonic wiggles (k ~ 0.1 Mpc -1 ) will be in this 
large scale regime. This is a region in k-space that con¬ 
tains much cosmological information and is also on scales 
where interferometers will be most sensitive to the p de¬ 
composition of the signal (Fig. 1101) . 

7. DISCUSSION 

In this paper, we have used the FZH04 analytical 
model of reionization to calculate the power spectrum 
of 21 cm brightness temperature fluctuations Pat(}+), 
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TABLE 4 

Errors on cosmological parameters estimates from 21 cm observations when only P„ 4 and P,,b 

/X f-L 

ARE NOT CONTAMINATED BY THE BUBBLES. a 



T 

flA 

flm h 

fib h ? 

n s 

5 h x 10 s 

Cts 


Hi 

"bT 

II 

00 


0.1 

0.7 

0.14 

0.022 

1.0 

3.9 

0.0 

0.0 

0.8 

CCMB 

0.060 

0.084 

0.017 

0.0014 

0.07 

0.3 

0.039 

0.12 

- 

CCMB + MWA5000 

0.058 

0.066 

0.011 

0.0012 

0.05 

0.2 

0.033 

0.05 

0.3 

CCMB + SKA 

0.057 

0.062 

0.011 

0.0012 

0.04 

0.2 

0.028 

0.04 

0.3 

Planck 

0.005 

0.011 

0.0023 

0.00017 

0.0047 

0.03 

0.007 

0.010 

- 

Planck + MWA5000 

0.005 

0.011 

0.0022 

0.00017 

0.0047 

0.03 

0.007 

0.010 

0.07 

Planck + SKA 

0.005 

0.011 

0.0022 

0.00017 

0.0047 

0.03 

0.007 

0.010 

0.07 

Planck + MWA50K 

0.005 

0.010 

0.0020 

0.00017 

0.0047 

0.03 

0.007 

0.009 

0.06 


a See table note (a) in Table|2]for the specifications used in these calculations. Observations with SKA and 
MWA50K are for 400 hours on ten places in the sky. 


extending this calculation beyond calculations done in 
FZH04 by including redshift-space distortions. When 
i, < 0.5 or on smaller scales than the effective bubble 
size, these distortions are quite important and give the 
power spectrum a substantial anisotropy between modes 
parallel and perpendicular to the LOS. These distortions 
not only increase the signal, but may allow us to sep¬ 
arate P xx , P x s and Pss, facilitating the measurement of 
the size and bias of the HII bubbles and perhaps the spec¬ 
trum of density fluctuations in the Universe. We show 
that higher order terms may complicate the separation 
of Pat when the ionized fraction is significant. 

To quantify the detectability of Pat(I^), and in par¬ 
ticular Pgs, we make realistic sensitivity estimates for 
LOFAR, MWA and SKA. The most important parame¬ 
ter for these interferometers is the collecting area. But, 
this is not to say that the other factors that go into the 
design a re unim po rtan t. We agree with the conclusion of 
iBowman et al.l 112005all that, everything else being equal, 
arrays with denser cores will be more sensitive to Pat{ k). 
This is because modes along the LOS can be detected by 
even the shortest baselines, and arrays with cores have 
more of these shorter baselines. The antenna size can 
also have a similar effect: arrays that have large anten¬ 
nae cannot pack them as closely together as arrays with 
small antennae. As a result, they will not sample the 
shorter baselines as well. Smaller antennae also provide 
a larger FOV, which will aid statistical detections of the 
signal. Because the current design for LOFAR does not 
include the shorter baselines that the design for MWA 
has and because the design for LOFAR results in a much 
smaller FOV, we find that, despite differences in collect¬ 
ing area, LOFAR and MWA will be comparably sensitive 
to Pat(^) at most redshifts. This is not to say that this 
will be the case when these instruments are actually de¬ 
ployed. Since none of the discussed 21 cm arrays have 
begun construction, their designs can still be optimized. 

Even with an optimally constructed radio interferome¬ 
ter, the removal of foregrounds that are 10 4 times larger 
than the 21 cm fluctuations will be a serious challenge. In 
this paper, we find that foregrounds will contaminate the 
signal on scales greater than the depth of the slice used 
to construct the 21 cm power spectrum. At most scales 
smaller than this, we are optimistic that foregrounds can 
be cleaned below the signal. On such scales, we show that 
fitting a quadratic or cubic polynomial to the observed 


visibilities in the frequency direction has little difficulty 
cleaning a realistic model for 
pea r to be the case, as was clc 
an d lGnedin fc Shaver! 120041) 
taminate all angular modes beyond repair. 

Applying our calculation for the detector noise and 
foreground power spectrum, we find that MWA and 
LOFAR will not be sensitive to the P ^4 component of 
PAT(k). This component is particularly interesting be¬ 
cause it traces the linear-theory density power spectrum. 
However, these interferometers will be sensitive to P^o, 
which probably will tell us more about the astrophysics 
of reionization than about cosmology, except perhaps for 
very small ionization fractions. MWA5000 and SKA will 
be moderately sensitive to P^i and P^a , but not sensitive 
enough to provide competitive constraints on cosmolog¬ 
ical parameters. We find that only if there exists a time 
when density fluctuations dominate Pat, will upcom¬ 
ing probes of 21 cm emission be able to place competi¬ 
tive constraints on cosmological parameters. In addition, 
planned 21 cm interferometers will not be very sensitive 
to the signal for 2 > 12. The is primarily because detec¬ 
tor noise fluctuations are proportional to T s k y (.z), which 
scales as (1 + z) 26 . 

If there is a period where density fluctuations dominate 
Pat, a 2 yr observation with MWA5000 plus Planck can 
give the constraints Sfl w = 0.0017 (a 1.7 times smaller 
uncertainty than from Planck alone), Sw = 0.06 (1.5 
times), 6Q m h 2 = 0.0009 (2.5 times), Sflbh 2 = 0.00012 
(1.5 times), Sn s = 0.0033 (1.4 times), Sa s = 0.003 (2.7 
times), 8£l v = 0.003 (3 times) and Sxh = 0.03. SKA 
plus Planck yield similar constraints and MWA50K can 
do even better. However, if r = 0.17, as suggested by 
WMAP, and reionization began at z ~ 20, observations 
of the signal at scales much larger than the effective bub¬ 
ble size may be the most promising direct method to 
probe cosmology m . 

Observations must overcome many additional chal¬ 
lenges beyond those that have been discussed in this 
paper. Issues that we have not addressed include con¬ 
tamination by radio recombination lines, terrestrial ra¬ 
dio interference, the residuals left from wave front correc¬ 
tions for a turbulent atmosphere, and the enormous data 
analysis pipeline needed to analyze potentially larger 
data sets than those from current experiments. The 
21 cm signal will also be affected by gravitational lens- 


foregroun ds. I t doesji ot^ap- 
, that foregrounds will con- 
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ing b y intervening material liZahn fc Zaldarriagal l2005t 
iMandel fc Zaldarriagall2005lh If taken into account, this 
effect can furt her imp r ove cosm ological parameter esti¬ 
mates liZahn fc Zaldarriagall2005th 
Cosmic variance sets a limit on how well we can 
constrain cosmological models with the CMB. Because 
21 cm emission can be observed as a function of redshift, 
this signal allows us to measure many more independent 
modes than is possible with the CMB. We have seen 
that cosmological parameters are extractable from 21 cm 
emission. In an ideal case in which reionization begins at 
relatively low redshifts, upcoming interferometers may 
be able to compete with future CMB experiments such 
as Planck. If reionization begins at higher redshifts or if 


the spin temperature fluctuations are important, a more 
sensitive interferometer will be required than those that 
are currently planned to be able to compete with CMB 
parameter constraints. Regardless of how reionization 
actually proceeded, high-redshift 21 cm emission has the 
potential to become a valuable probe of cosmology. 
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71019, AST 02-06299, and AST 03-07690, and NASA 
ATP grants NAG5-12140, NAG5-13292, and NAG5- 
13381. 


APPENDIX 

THE EFFECT OF EVOLUTION 

Observations must have a large enough bandwidth to provide adequate signal-to-noise, but the larger the bandwidth, 
the more the signal will evolve along the LOS direction of the 21 cm map. During reionization this evolution includes 
( 1 ) density inhomogeneities growing with time and ( 2 ) the average 21 cm brightness temperature declining as the 
Universe expands and the bubbles grow and occupy progressively more space. The latter should dominate over the 
former owing to the relatively short timescales during which the ionization fraction changes by order unity. Since the 
Universe is evolving across modes along the LOS and is not evolving for those transverse to the LOS, this evolution 
will break the spherical symmetry of the signal and potentially make it more difficult to observe P^o , 2 , P ^4 and 

P^e . In this section, we attempt to quantify the potential size of this effect for an idealized case. 

Very simple models for reionization set Pat = 6 2 P$$ where b is the effective bias. This parameterization is certainly 
not correct, and the bias will have a scale dependence for scales around the size of the bubbles. This model for the 
power spectrum is reasonable at the beginning of reionization, when density fluctuations dominate, and on much larger 
scales than the bubbles . Fortunately, these are the regimes in which the cosmological information is most readily 
extractable from the signal (see eh and 

With the assumption that Pat = b 2 P 55 , we can write an expression for the power spectrum in a region of width 
2 Ar 

PAr{k) = f d 3 rd 3 eF(r ll ) 6 0 (r)F(r ll +e l \) 6 0 (r + e)exp(-ik-e)), (Al) 

Al/ J AV 

where AV is specified in the LOS direction by the conditions |v|| — f|j| < Ar and |r|| + e|| — f||| < Ar and r|| is the 
projection of r along the LOS and is assumed to be much larger than Ar in the angular directions.. To constants of 
order unity, F = b(z ) (1 + z) -1 / 2 [the factor of (1 + z) -1 / 2 owes to the evolution of Xb as well as the growth factor], 
and So are the linear density fluctuations at z = 0 [Sl(z) = G(z ) <5o, where G is the growth factor]. Here the subscript 
|| indicates the LOS direction. The () indicates an ensemble average of maps, and (So(r)So(r + e)) = C<5o<5o( e ) i n which 
£<5 0 <5 o is the linear density field corr elation function. The Fourier transform of is Pg 0 s 0 - By the Fourier transform 
convolution identity, equation m is equivalent to 

P AT (k) = jdk'^PsMk-k'), (A2) 

where 

W(k\\) = [ dr \\ de \\ F ( r \\) F ( r \\ + e ||) ex P(-*fc|| £||)- (A3) 

Jav 1{ 

Prior to reionization, 6 = 1 and W{k) is almost identical to the window function for a top-hat in real space. Since 
this window function peaks when its argument is near zero, for large k , then PAr(k) ~ Pg 0 s 0 (k) J dk'W(k'). This 
is exactly what we should expect; in this limit, we are effectively averaging over modes located at different redshifts, 
such that only the amplitude of the observed power spectrum is affected by evolution. The difference between this 
window function for a universe in which the ionized fraction is and is not changing will indicate the degree to which 
the spherical symmetry is affected by the evolving ionized fraction. For example, if the window function becomes less 
peaked, then the spherical symmetry of more modes is affected by the evolving ionized fraction. In cases in which the 
window function is unaffected, this means that evolution only affects the normalization of the power spectrum through 
a factor of J dk' W(k'). While the normalization is also interesting, it is degenerate with the parameters xh , 61 (the 
linear bi as of t he bubbles) and erg and therefore is less important. 

Figure EH plots the window function at z = 10 for maps with comoving depth 110 and 230 Mpc (Az = 0.5 and 
1) for a neutral universe (solid curve) and two cases in which 6(z) oc [m(z — 10) + 1] for m = 0.5 and 1. This is 
a conservative choice; such a quick change of b is much larger than models predict at the beginning of reionization. 
Figure EH demonstrates that the window function is virtually unaltered for three of the four cases. Only in the most 
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Fig. A12.— Window function W(k), defined in equation <A3l . for maps of comoving depths 110 Mpc ( outer curves) and 230 Mpc ( inner 
curves) at z = 10 . We plot the window function for a universe prior to reionization b(z) = constant ( solid curves) and for when the 
universe is starting to be ionized (quickly) with b(z) oc \rn{z — 10) + 1] for m = 0.5 and 1 ( dot-dashed and dashed curves, respectively). 


extreme case, where m(z — 10) + 1 varies from 0 to 1 within the map, is the window function significantly altered. 
Therefore, the evolution of the signal for this simple model weakly breaks the spherical symmetry of the signal. 


ANGULAR AVERAGED SENSITIVITY 


The spherically averaged signal is obtained by summing up all pixels in a shell with the same \k\. If A k = e/c, it 
follows from equation itfl) that the error on .PAT(fc) from a measurement of all pixels in a shell with constant k is 
given by 


SPAr(k) = 


E 


( 5-Pat (fc, 0 ) 


- 1/2 


/*arcsin[min( , 1)] 

k 3 / d,6 sin(0) 

a arccos[min(^,l)] 


1 - 1/2 


D = 


(2tt) 2 A e 
A 2 x 2 ye 


E = 


(p P21cm(^) + „( fes f n (0)) 
2tt Pxhjx 3 r s 2 ys 


(Bl) 

(B2) 


where fc* corresponds to the longest wavevector perpendicular to the LOS probed by the array. In equation m, the 
lower bound of the integrand reflects the sharp cutoff in the number of pixels imposed for wavelengths that do not 
fit in the width of the box, y. Some interferometers, such as LOFAR, will be able to observe N po [ nt separate fields of 
view simultaneously. To include this effect, one can simply divide equation m by ^/iVp 0 i n t. 

To gain intuition into the scalings of equation m , let us take a top-hat distribution of baseli nes with density p( A) 
in the limit in which the detector noi se d ominates such that we can set P 2 icm = 0 in equation m■ In this case, we 
can evaluate the integral in equation EQ>, which yields 


SPAr(k)=E ^min[|^,l] 

\“ 1/2 

1 p- 1 k ~ 3/2 

k < fc* 

(B3) 

= Ep- 1 k~ 3/2 { 1 

\Jk 2 - k 2 1/2 
k 1 

k > k*. 

(B4) 


When k fc*, <5 Pat(^) oc k 0 5 . 
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